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1  Introduction 


The  accuracy  of  a  design  and  a  cost  estimate  depend  upon  the  stage  at  which 
they  are  performed.  This  accuracy  ranges  from  5C  at  the  definitive  stage, 
when  almost  all  the  engineering  data  and  design  details  about  the  project 
are  available,  to  more  than  35c~  at  the  feasibility  stage,  when  there  is  less 
information  available  about  the  project  9  The  objective  here  is  to  develop 
a  model  that  produces  a  reliable  basis  fur  estimating  at  the  feasibility  stage 
This  stage  is  very  important  for  both  the  owner  and  the  designer.  It  is 
necessary  to  have  a  rebable  basis  for  design  and  an  accurate  cost  estimate  at 
this  stage  in  order  to  decide  whether  to  proceed  with  the  project  or  net. 

To  accomplish  this,  a  review  of  the  types  cf  estimation  methods  available 
is  first  presented.  Two  methods,  using  regression,  that  deal  with  the  estimate 
at  the  predesign  stages,  are  also  reviewed.  In  such  estimates  there  are  many 
subjective  factors  that  can  be  handled  by  using  the  fuzzy  sets  concept  A 
mode!  for  analyzing  design  reliability  based  on  fuzzy  sets  theory  is  proposed 
Numerical  examples,  including  sensitivity  analysis,  are  also  presented 


2  Types  of  Estimation 

Cost  estimates  range  from  a  feasibility  study  to  complete  design  of  project 
The  estimate  and  type  of  design  data  desired  may  differ  radically  from  one 
project  to  another  In  this  section,  a  renew  of  the  available  methods  at  all 
stages  will  be  given  and  the  pitfalls  ■:>!  the  estimates  will  be  outlined  A 
detailed  review  of  the  preliminary  estimate  will  fc.liow 

2.1  Estimates  and  their  accuracies 

As  mentioned  previously,  several  cost  estimates  are  usually  performed  at 
different  stages  of  any  project.  The  accuracy  of  an  estimate  depends  upon 
the  reliability  of  the  design  which  is  performed  Usually,  the  stages  and  their 
corresponding  accuracies  are  as  follows 

1.  Definitive  =s  accuracy  ro'c 

2.  Appropriation  =>  accuracy  nK1  -  1 5°T 
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3  Feasibility  =>  accuracy  >  —  35  u 

It  is  very  important  to  note  that  the  estimate  in  all  stages  is  an  approxi¬ 
mation  based  on  judgment  and  expeiience.  Some  of  the  reasons  that  lead  to 
inaccurate  estimates  are  pointed  out  in  Section  2.2. 

2.2  Pitfalls  of  the  Estimation 

The  pitfalls  of  the  estimation  may  vary  from  one  project  to  another  The 
following  pitfalls  may  be  common  to  many  projects  9  . 

1.  Misinterpretation  of  the  statement  of  work. 

2.  Omissions  or  improperly  defined  scope. 

3  Poorly  defined  or  overly  optimistic  schedules. 

4.  Inadequate  project  design  data  and  engineering  properties. 

5.  Inaccurate  work  breakdown  structure. 

6.  Inappropriate  analytical  techniques  or  modeling. 

7.  Applying  improper  skill  levels  to  tasks. 

8.  Failure  to  account  for  risks,  including: 

•  safe  and  reliable  design  concepts 

•  site  geology 

•  labor  productivity 

•  unproven  technology 

•  regulatory,  environmental  factors 

•  equipment  productivity 

•  changed  conditions 

•  permit  approval 

•  weather 

•  accidents 
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•  quality  control  during  design  and  construction 
9  Failure  to  understand  or  account  for  cost  escalation  and  inflation. 

10.  Failure  to  use  correct  estimating  techniques. 

11  Failure  to  use  forward  pricing  rates  for  overhead,  general  and  adminis¬ 
trative,  and  indirect  costs. 


2.3  Preliminary  Estimate 

A  preliminary  estimate  is  generally  presented  before  all  the  information  re¬ 
quired  for  the  project  is  available.  Such  estimates  are  performed  by  the  owner 
or  his  consultant  during  a  feasibility  study,  by  the  designer  to  evaluate  pos¬ 
sible  design  alternatives,  and  by  contractors  for  bidding  and  budgeting.  A 
cost  estimate  at  the  feasibility  study  stage  cannot  be  as  precise  since  it  is 
based  on  an  approximate  design.  However,  an  accurate  method  is  required. 
To  develop  such  a  method,  the  identification  of  available  methods,  allocation 
of  risk  involved  as  well  as  when  it  is  likely  to  occur,  and  identification  of 
design  and  cost  variables  responsible  for  large  variations  in  the  project  cost 
are  required. 

2.3.1  Types  of  Preliminary  Estimate 

Most  of  the  existing  conceptual  and  preliminary  est;  ating  methods  fail  into 
one  or  more  of  the  following  categories  3  : 

1.  Time  referenced  indices. 

2.  Capacity  factors. 

3.  Component  ratios. 

4.  Parameters. 

These  methods  will  be  explained  in  detail  in  sections  2  3.2  to  2.3.5. 
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2.3.2  Ind  ices 


Indices  show  changes  ir.  estimated  quantities  over  time  These  indices  are 
divided  into  two  types.  The  first  type  updates  the  data  which  serve  as 
inputs  to  the  project.  The  second  one  is  based  on  the  completed  construction 
project,  which  is  based  on  the  output  of  the  project. 

2.3.3  Capacity  Factors 

Capacity  factors  apply  to  changes  in  size,  scope,  or  capacity  of  projects  of 
similar  type,  as  opposed  to  indices  which  focus  on  changes  over  time.  Time- 
independent  factors  are  expressed  in  the  following  equation  3  : 

C2  =  (7,(52)*  (!) 

Vi 

where  C 2  is  the  estimated  quantity  of  the  new  project  of  capacity  02.  Cj  is 
known  quantity  of  facility  of  capacity  Qx.  and  x  is  the  capacity  factor  for 
this  type  of  work. 

Capacity  factors  and  indices  can  be  combined  to  take  into  account  changes 
in  both  time  and  capacity.  The  modified  time-dependent  formula  is  3': 

7,  0, 

Cl  =  (2) 

■<i  0i 

where  C%  is  the  estimated  quantity  of  the  new  project  and  /j  and  12  are  the 
relative  indices  for  the  time  associated  with  known  and  proposed  facilities, 
respectively. 

2.3.4  Component  Ratios 

The  component  ratio  method  is  based  on  information  concerning  the  major 
equipment  to  be  installed  in  the  project.  This  can  be  found  by  multiply¬ 
ing  the  equipment  characteristic  or  property  by  an  empirically  documented 
factor. 

2.3.5  Parameters 

The  parameter  approach  relates  all  design  aspects  of  a  project  to  just  a  few 
physical  measures  of  parameters  that  reflect  the  size  or  scope  of  that  project 
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For  example,  the  gross  enclosed  floor  area  would  be  a  typical  parameter  for 
a  structure  such  as  a  warehouse. 


3  Regression  Methods 

3.1  Kouskoulas  Method[20] 

In  using  this  method,  the  building  cost  function  is  defined  as  a  function  of  six 
variables,  namely,  locality  index,  price  index,  building  type,  building  height, 
building  quality  and  building  technology.  The  data  used  in  this  method  are 
also  used  in  the  proposed  model,  therefore,  this  method  will  be  exolained 
later  on  in  Section  4. 


3.2  Karshenas  Method[l7] 

In  the  Karshenas  (1984 )  17  model,  the  equation  selected  is  a  power  equation. 
The  cost  in  this  model  is  a  function  of  height  and  floor  area  only,  however, 
these  are  not  the  only  variables  that  affect  the  cost. 

3.2.1  Data  Collection 

Engineering  News  Record(ENR)  data  since  1968  are  used  in  Karshenas' 
study.  Table  1  shows  the  data  for  office  building  construction  costs  Since  the 
year  of  construction  and  the  location  of  each  building  varies,  their  cost  cannot  be 
directly  compared  to  what  was  reported  by  ENR  In  order  to  make  the  building 
costs  comparable,  all  costs  were  converted  to  March  1982.  New  York  City  costs,  by 
using  the  cost  index  for  steel-framed  buildings,  published  by  Marshall  and  Swift 
Publishing  Company  17  .  The  converted  costs  are  shown  in  the  second  column  of 
Table  2. 

3.2.2  Cost  Function 

The  data  in  Table  1  may  be  used  to  derive  a  cost  function  of  the  general  form 

C  =  f  { A,  H)  (3) 

in  which  .4  and  H  are  the  typical  floor  area  and  height  of  the  building,  respectivelv 
One  method  for  investigating  the  shape  of  this  cost  function  is  n  pb  t  contours 
of  constant  C  using  .4  and  H  as  abscissa  and  ordinate,  respectively  Figure  1  shows 
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an  example  of  such  contours  fitted  "by  eve".  The  least  squares  method  was  used 
to  select  the  type  of  function  and  power  function  was  chosen.  The  general  form  of 
the  power  function  is  1 7  : 

f{A.H)=A3H "  (4) 

in  which  3  and  q  are  constants.  Therefore,  Eq.  (4)  would  be 

C  =  aA5H '  (5) 

in  which  a  is  a  constant 

Taking  the  natural  logarithm  on  both  sides.  Eq.(5)  may  be  equivalently  ex¬ 
pressed  as 

lnC  =  LnQ-51n.4-7ln/f  (6) 

Using  the  data  given  in  Table  1  and  the  least  squares  method,  the  parameters 
of  the  preceding  equation  were  estimated.  Eq  (6)  with  the  estimated  parameters 
is 

In  C  =  -0.0235  -  1.1045  In. 4  -  1.1268  In#  (7) 

or 

q  _  ^1.1045  ^  1 .1268^-0. 0235 

in  which  C  is  the  predicted  cost  of  a  budding  with  the  typical  floor  area  A  (sq. 
ft),  and  height  of  H  (ft).  Examples  of  the  contours  of  Eq.(8)  superimposed  on 
the  scatter  plot  of  the  observed  budding  costs  are  shown  in  Figure  2  The  costs 
predicted  by  Eq.(8)  are  tabulated  in  Column  3  of  Table  2 


3.2.3  Examination  of  the  Formula 

The  fodowing  questions  are  considered  in  this  section: 

1.  Whether  the  data  are  adequately  described  by  the  regression  equation'’ 

2.  How  accurately  Eq.(8)  predicts  the  cost  of  a  budding"’ 

The  ability  of  the  model  to  explain  the  variations  of  the  cost  data  can  be 
measured  by  investigating  the  coefficient  of  determination,  defined  as 


where 


St  =  Dc>  -  n2 

i=i 

and 

n 

Sr  —  ^  J  Ci  Uo  ~  i|  —  .  Uml  mi)" 

1  =  1 


6 


Figure  2:  Scattergrara  of  Observed  Costs  and  Contours  of  Estimated  Regres¬ 
sion  Plane  (in  millions  of  dollars)  [17] 
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where  C  is  the  mean  value  of  the  cost  and  v  is  the  number  of  observed  data 

Lf  Eq  (8)  perfectly  predicts  the  cost  of  the  building,  (i.e.,Sr  =  0)  .  then  r  =  1. 
On  the  other  hand,  if  the  proposed  model  does  not  explain  the  variations  in  the 
building  cost,  r  =  0. 

4  Feasibility  Estimate 

The  design  of  any  project  is  a  function  of  many  variables.  An  important  problem 
to  consider  is  the  selection  of  a  set  of  independent  variables  that  defines  a  project 
design  and  its  cost.  Another  problem  is  the  selection  of  the  cost  function  in  terms 
of  the  design  variables.  The  criterion  for  selecting  the  variables  is  availability  of 
data  concerning  such  variables  from  previously  completed  projects.  Regarding 
data  availability,  the  cost  function  may  be  defined  as: 

C  =  f(Vl, . ,Vm)  (10) 

in  which  C  is  a  cost  measure  of  a  project  and  1',  are  independent  design  variables. 
The  design  variables  may  include: 

1.  Location,  including  site  geology. 

2.  Quality  of  work. 

3.  Labor  productivity. 

4.  Material. 

5.  Equipment  productivity. 

6.  Technology. 

7.  Weather. 

8.  Project  type. 

9.  Structure  height. 

4.1  Cost  Function 

The  selection  of  independent  design  variables  and  the  construction  of  the  cost  func¬ 
tion  depend  on  the  availability  of  data.  The  design  variables  used  by  Kouskoulas 
(20;  are  the  data  for  a  building.  His  independent  design  variables  include,  building 
locality  index,  price  index,  building  type  index,  building  height,  building  quality 
index  and  technology  index. 
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4.2  Measurement  of  Design  Variables 

In  preliminary  cost  estimating,  one  may  identify  a  building  by  its  location,  time  of 
realization,  function  or  type,  height,  quality  and  technology.  Consider  the  variable 
l'j  =  L.  identifying  the  differences  in  construction  costs  as  a  consequence  of  dif¬ 
ferences  in  the  life  style  between  different  cities.  Tabie  3  presents  typical  locality 
indices  for  ten  cities.  The  price  index  variable  1  \  —  P  is  time- dependent.  The 
problem,  therefore,  is  that  the  price  index  has  to  be  projected  into  the  future  or 
to  the  time  when  the  construction  takes  place  in  order  to  construct  the  function. 

P  =  g(t)  (ii) 

Utilizing  historical  data,  Table  4  provides  a  sample  time  history  of  average 
construction  costs  in  dollars  per  cubic  foot  for  various  apartment  buildings. 

The  third  variable,  U3  =  T,  specifies  the  types  of  buildings.  A  good  measure 
for  this  variable  is  provided  by  the  average  cost  per  cubic  foot  for  different  types  of 
buildings.  Table  5  shows  a  range  of  classes  of  buildings,  with  their  corresponding 
relative  cost  values,  as  provided  by  [20). 

The  height  index.  U4  =  H .  indicates  the  number  of  stories. 

The  quality  variable.  U5  =  Q ,  stands  for  what  it  specifies  It  is  a  measure  of: 

1.  The  quality  of  workmanship  and  materials  used  in  the  construction  process 

2.  The  building  use. 

3  The  design  effort. 

4.  The  material  type  and  quality  used  in  various  building  components. 

To  define  the  measure  for  such  a  variable,  let  C,,i  =  1.  ...k,  denote  the  average  cost 
portions  of  total  building  cost  distributed  among  k  building  components.  Assign 
an  integer  index,  I  ,  from  1  to  4  (corresponding  to  fair,  average,  good  and  excellent) 
to  each  component,  and  compute  the  values  of  the  variable,  Vs,  from[20  : 

li.jEw  (i2) 

K  .=i 

A  more  subjective  way  might  be  to  identify  the  components  of  the  building, 
assigning  each  a  value  of  1  to  4  and,  on  this  basis,  pass  a  judgment  on  the  quality 
of  the  total  building  by  assigning  it  one  of  the  integer  values  1  to  4.  It  must  be 
realized  that  the  range  of  integer  value.  1  to  4.  is  arbitrarily  selected  and  is  only 
good  as  long  as  it  can  provide  a  reliable  estimation  function. 
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Table  6  identifies  general  building  components  and,  on  the  basis  of  their  qual- 
itat've  description,  rates  them  accordingly,  from  fair  to  excellent. 

The  technology  index  variable.  I’6  =  TE,  takes  into  account  the  extra  cost 
that  is  incurred  for  special  types  of  buildings  or  the  labor  and  material  savings 
that  result  from  the  use  of  new  techniques  in  the  construction  process.  Some  data 
regarding  this  variable  are  given  in  Table  7. 


4.3  Derivation  of  Cost  Function 

The  cost  function.  C  =  may  take  many  possible  forms.  The  one 

chosen  here  is  the  linear  expression: 


C  -  a0  -  Ul  1 1  -  <22!  2  — . —  Uml  m  ~  e  (13) 

where  afc(Jfc  =  0 ,....,m)  are  constants  to  be  determined  from  the  available  data 
and  (e  =  C,  -  C,)  is  the  error  associated  with  the  selected  function.  Once  the 
coefficients  a ^  are  known,  C  can  be  computed  for  any  new  project  identified  by  a 
set  -f  -r.iues  Vk  .  The  coefficients  a*  (Jt  =  0,  l,...m)  are  found  by  multiple  linear 
regression  analysis  using  the  data  in  Table  8  . 
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Table  1:  Office  Building  Design  Data  and  Costs  17 


Number 

Location 

Time  of 

construction 

Number 

of 

floors 

Building 

height 

(feet) 

Typical  floor 
area 

(square  feet) 

Total  j 

Cost 
(dollars) 

1 

Lexington, Mass 

Nov  77/Jan.  79 

3 

36 

27,000 

3,133,100 

2 

Southfield, Mich 

Apr  76/June  78 

15 

187  5 

17,690 

11,645,000 

3 

Pocatello,  Id  alio 

May  76/Sept  77 

3 

36  6 

26,690 

2,969,800 

4 

Dallas, Tex 

Apr  76/May  77 

21 

262  5 

17,000 

12,958,000 

5 

Glen dole, Calif 

Dec.75/'Nov  76 

6 

72 

15,800 

2781,480 

6 

Seattle,  Wash. 

Feb  74 /Dec  76 

36 

468 

22,200 

36,470,000 

7 

Scottsdale,  Am. 

Dec  74/May  76 

2 

28 

140,332 

1 2,729,90u 

8 

Knoxville, Tenn. 

Aug  74/Apr  75 

2 

25  3 

9,986 

446,500 

9 

Troy, Mich. 

Aug  73/Oct  75 

26 

330 

19,400 

16,822,00 

10 

Birmingham, Ala. 

Oct  74/Jan. 76 

18 

216 

12,616 

6,104,140 

11 

Franklin,  III 

Mar  74/dec,74 

5 

62  6 

8000 

1,396,200 

12 

Beverly  Hills, Calif. 

Nov  73/July  75 

8 

105  6 

5,500 

1,204,100 

13 

Houston, Tex. 

July  73/Jan. 75 

13 

175.5 

29,920 

10,408,000 

14 

Chicago,  Ill. 

Dec  73/Dec.  74 

2 

28 

35,280 

1,951,175 

15 

Detroit,  Mich. 

Aug  71/Apr  73 

4 

48 

17,700 

1,731,800 

16 

Warren, Mich 

June  72/Oct. 73 

11 

137  5 

15,000 

4,435,000 

17 

Wellesley,  Mass. 

Dec.  69/Sept  70 

4 

48 

18,800 

1,763,000 

18 

Central, NJ. 

Nov  70/Feb  72 

12 

153 

30,134 

1  1,129,000 

19 

San  Francisco, Calif 

Oct. 66/May  68 

33 

429 

17,212 

14,455,000 

20 

New  York,  N  Y. 

Oct  61  /Nov  63 

42 

483 

18,893 

16,820,900 

21 

Cleveland,  Ohio 

Feb.  63/Nov  64 

41 

533 

21,600 

20,1 16.000 

22 

Columbus,  Ohio 

Dec.  63/Feb  65 

26 

338 

16,000 

8.683,000 

23 

Pittsburgh,  Pa 

Apr  66/ Apr  68 

9 

126 

16,833 

3,871,000 

24 

Houston, Tex. 

Sept  65/'Aug  66 

1 

50 

10,500 

656,100 

12 


Table  2: 

Predicted  Cost 

and  Residu 

als 

Num¬ 

ber 

Adjusted  total 
cost.c  . 
(dollars') 

Predicted  total 

cost  ,c. 
(dollars') 

Residuals, 

e  =  c  -  c 

Residual. as  a 
percentage.'e  c) 

1 

4.542.000 

4.377.454 

164.545 

0  036 

2 

17.271.000 

17.619,844 

-348.644 

-0.020 

3 

4.699,000 

4.441.512 

257.487 

0.054 

4 

25.254,000 

24.636.604 

i  617.396 

0.024 

5 

4.506,000 

5.289.436 

j  -783.436 

-0.170 

6 

70.022.000 

63.469.742 

6.552.258 

0.093 

t 

21.050.000 

20.362.490 

687.510 

0.032 

8 

928.720 

980.603 

-51.883 

-0.055 

9 

31.942.800 

36.890.823 

-4.946.023 

-0150 

10 

12.704.000 

14.226.597 

-1.522.597 

-0.120 

11 

2.619.110 

2.126.700 

492.410 

0.190 

12 

2.175.000 

2.538.924 

-363.924 

-0.160 

13 

23.722.000 

29.222.684 

5.500.684 

-0.230 

14 

3.902.350 

4.431.366 

-529.016 

-0.130 

15 

3.300.000 

3.797.100 

-497.100 

-0.150 

16 

9.275.000 

10.353.680 

-1.078.680 

-0.110 

17 

4.531.000 

4.058.570 

472.429 

0.100 

18 

26.932.000 

25.234.518 

1.697.462 

0.06 

19 

42.931,000 

43.442.324 

-511.324 

-0.012 

20 

63.346,000 

55.034.241 

8,311.759 

0.130 

21  1 

69.946.000 

71.296.732 

-1.350.732 

-0.019 

22  i 

31.817,000 

30,634.697 

1.182.303 

0.037 

23  j 

11,818.000 

10,657,469 

1,160.531 

0.098 

24 

2,646.100 

2,233,291 

412.809 

0.15 
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Table  3:  Typical  Locality  Indices  for  10  Cities  20 


City 

Index. L 

i  Boston, Mass. 

i  1.03 

Buffalo. NY. 

1  1.10 

i  Dallas. Tex. 

0.87 

Dayton. Ohio 

i  1.05 

'  Detroit. Mich. 

1.13 

1  Erie. Pa. 

1.02  ; 

Houston. Tex. 

0.90 

Louisville.  Ky. 

1  0.95  ■ 

New  York  .N  Y. 

1.16 

Omaha. Neb 

0.92 

Table  4:  Sample  Time  History  of  Average  Construction  Costs,  in  Dollars  per 
Cubic  Foot, for  Various  Apartment  Buildings  20 


Year 

t 

Index. P 

1  1963 

o 

1.66  ; 

1964 

1 

1.71  i 

1965 

2 

1.76 

1966 

3 

1.8  ! 

1967 

4 

1.93  j 

1968 

5 

2.09 

1969 

6 

2.30 

1970 

7 

2.49 

i  1971 

8 

2  76 

|  1972 

9 

2.95 

14 


Table  5:  Range  of  Classes  of  Buildings  with  Corresponding  Relative  Cost 
Values  20 


T  vpe 

Index.! 

Apartment 

2.97 

Hospitals 

3.08 

Schools 

2.59 

Hotels 

3.08 

Office  building!  fireproof  I 

2.95 

Office  building! not  fireproof! 

1  83 

Stores 

2  43 

Garages 

1.99 

Factories 

1.2 

Foundries 

1  49 

15 


Table  6:  Quality  Index  0  20 


Component 

Fair 

Average 

Good 

Very  Good 

Use 

Multitenancy 

1 

Mixed,  single 
tenant .  and 
Multitenancy 

Single  tenant 

i 

Single  tenant 
with  custom 

requirements 

Design 

Minimum 
design  load 

Average 
design  load 

Above  average 
design  load  1 

Many  extra 
design  load 

Exterior  Wall 

Masonry 

Glass  or 

masonrv 

1 

Glass. curtain 

wall, precast 

1  1 
concrete 

Monumental 

(marble) 

Plumbing  i 

! 

Below  Average 
quality 

Average 

quality 

Above  average 
quality 

Above  average 
qualit  v 

Flooring  1 

1 

! 

1 

- 1 

Resilient, 

ceramics 

Resilient. 

ceramics 

and 

terrazzo 

Vinyl. ceramic 

terrazzo 

Rug.  terrazzo 
marble 

Electrical 

1 

1 

1 

Fluorescent 
light, poor 
quality  ceiling 

Fluorescent 
light  .average 
quality, 

suspended  ceiling 

Fluorescent 
light,  above 
average 
quality  ceiling 

Fluorescent 
light  .excellent 
qualit  v 
ceiling 

i  Heating,  j 

Below  average 

Average 

Above  average 

Above  average 

!  Ventilating 
and  Air 
Conditioning 

quality 

quality 

i 

quality  | 

1 

quality 

|  Elevator 

i 

Minimum 

required 

i  Above  required 

i 

High 

speed 

High  speed 
deluxe 

16 


Table  7:  Data  Regarding  Magnitude  of  Vi  [20] 


Technology  i 

Index. TE  j 

Bank-monumental  work 

1.75  I 

Renovation  building 

0.50 

Special  school  building 

1  10  : 

|  Chemistry  laboratory  building 

1.45 

■  Telephone  building- blast  resistant 

1  60 

County  jails 

1.20 

Dental  school 

1.15 

Hospital  addition 

1.05 

Correctional  center 

1.20 

i  Home  for  aged 

1.10 

17 


The  basic  idea  of  the  least  squares  fitting  method  is  to  minimize  the  sum-of- 
squares  of  the  residuals.  Sr  .  as  in 

n  n 

Sr  =  Y1  f2  =  ~  a°  ~  ail  n  -  .  -  amrmi)2  (14) 

t=l  »-l 

Then  the  function  is  differentiated  with  respect  to  each  of  the  unknown  coefficients. 

^  =  -2  y  Vki{Ci  -  a0  -  aiVu  -  . -  amrm.)  (15) 

oak 

r0,  =  1.  k  =  0, . m  .  i=l. . n 

The  coefficients  yielding  the  mimmum  sum  of  squares  of  the  residuals  are 
obtained  by  setting  the  partial  derivatives  equal  to  zero  and  expressing  Eq.  (15) 
as  a  set  of  simultaneous  linear  equations. 

m  n 

DLr*a*)  =  Er«c<  (16) 

*=o  1  =  1 

Or.  in  matrix  form: 

n  £ri. 

£  yu  £  ri2, 

£Ui  £'2,'), 

£  1  mi  £  1  1 

where  C,  is  the  cost  corresponding  to  data  point  i  and  \\l  is  the  variable  corre¬ 
sponding  to  variable  k  and  data  point  i. 

Constants  were  computed  using  a  regression  subroutine,  which  is  included  in 
the  program  (as  listed  in  Appendix  A),  with  the  historical  data  in  Table  8  (from 

;2°i) 

According  to[20],  the  resulting  required  function  is: 

C  =  -84.97694  -  26. 67051’]  ~  11.23471 ’a  -  6.23531'3 
201221'4  -  5.03921 5  -  31.219091s  (18) 

The  coefficients  in  Eq.  (19)  are  slightly  different  from  Kouskoulas  s.  however, 
this  might  be  due  to  different  numerical  methods  used  in  the  analysis 
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Table  8:  Building  Historical  Data 


Cost  C,in 
dollars 
per 

square 

foot 

Locality 

Index 

V, 

Price 

Index 

V, 

type 

Index 

Vs 

Height 

Quality 

Index 

v5 

Technology- 

Index 

v6 

Variable 

combin¬ 

ations 

36.00 

0.90 

2.76 

2.95 

40 

2 

1.00 

1 

25.00 

0.87 

2.49 

2.95 

18 

1 

1.00 

2 

68.50 

1.02 

2.76 

2.95 

6 

4 

1.75 

3 

31.90 

1.03 

2.49 

2.97 

8 

2 

1.00 

4 

36.50 

1.10 

2.49 

2.59 

11 

3 

1.00 

5 

23.30 

1.13 

2.95 

2.95 

5 

2 

0.50 

6 

40.00 

0.95 

2.30 

3.08 

14 

4 

1.00 

7 

56.00 

1.13 

1.93 

2.95 

3 

4 

1.60 

8 

40.00 

1.05 

2.09 

3.08 

S 

4 

1.00 

9 

21.70 

1.13 

2.09 

1.99 

2 

1.00 

10 

42.00 

1.00 

2.76 

2.95 

< 

4 

1.00 

11 

45.81 

1.16 

2.49 

2.59 

1 

4 

1.10 

12 

62.00 

1.16 

2.95 

2.59 

7 

4 

1.45 

13 

85.00 

1.00 

2.95 

3.08 

6 

4 

2.25 

14 

47.50 

1.00 

2.49 

3.08 

7 

4 

1.15 

15 

34.30 

1.13 

1.93 

3.08 

3 

3 

1.10 

16 

37.00 

1.13 

2.76 

2.95 

24 

1 

1.00 

17 

31.90 

1.13 

2.30 

2.95 

10 

2 

1.00 

18 

40.00 

1.13 

2.30 

2.95 

22 

3 

1.00 

19 

49.50 

1.13 

2.95 

2.95 

27 

3 

1.00 

20 

36.20 

1.13 

2.09 

3.08 

10 

3 

1.00 

21 

24.00 

1.13 

2.76 

1.83 

1 

1 

1.00 

22 

38.80 

1.13 

2.09 

3.08 

1 

4 

1.05 

23 

20.00 

1.08 

1.93 

2.95 

4 

1 

1.00 

24 

1880 

1.13 

1.93 

2.59 

2 

1 

1.00 

25 

34.70 

1.13 

2.09 

2.95 

5 

3 

1.00 

26 

15.10 

1.13 

1.93 

1.83 

2 

1 

1.00 

27 

18.10 

1.13 

1.22 

2.59 

3 

2 

1.00 

28 

39.00 

1.13 

2.3 

3.08 

4 

2 

1.20 

29 

36.00 

1.13 

2.09 

3.08 

2 

2 

1.20 

30 

21.10 

1.07 

1.66 

2.59 

6 

2 

1.00 

31 

24.30 

1.07 

1.93 

2.5 

6 

2 

1.00 

32 

30.00 

1.13 

1.93 

2.59 

6 

3 

1.00 

33 

27.50 

1.13 

2.30 

3.08 

2 

1 

1.00 

34 

11.30 

1.00 

2.09 

1.49 

1 

1 

1  00 

35 

14.50 

1.02 

2.09 

1.20 

3 

2 

1.00 

36 

10.00 

1.05 

2.09 

1.20 

1 

1.00 

37 

14  75 

0.92 

2.3 

1.2 

1 

2 

1.00 

36 
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5  Fuzzy  Model 

Design  reliability  depends  upon  the  subjective  evaluation  and  judgments  of  experi¬ 
enced  designers.  Available  data  can  be  used  to  develop  the  functions  which  relate 
the  dependent  to  the  independent  design  variables.  With  the  regression  model, 
values  of  independent  design  variables  can  be  substituted,  for  example,  into  the 
cost  function  to  predict  the  cost  of  new  designs  This  can  be  done  if  we  have  pre¬ 
cise  information  about  these  variables  Unfortunately,  this  is  not  always  the  case 
In  many  cases,  precise  information  is  not  always  available.  The  only  thing  known 
about  some  variables  may  be  what  could  be  called  “verbal  descriptions.’"  such  as 
the  lateral  loading  may  be  ’’high,’"  the  location  is  “good.”  and  so  on.  Normally 
variables  like  these  are  omitted  from  the  model,  but  this  just  results  in  more  errors. 

The  objective  here  is  to  try  to  use  the  cost  function  developed  by  Kouskoulas 
Eq.  (18)  and,  rather  than  eliminating  independent  design  variables  which  are 
not  precisely  known,  use  fuzzy  set  theory  to  handle  these  types  of  variables  in  a 
regression  model. 


5.1  Fuzzy  Set  Theory 

A  conventional,  or  non-fuzzv  set  was  defined  by  the  use  of  an  indicator  function, 
for  which  any  element  of  the  sample  space  was  1  or  0.  depending  upon  whether  the 
element  was  a  member  of  the  set  or  not.  The  indicator  function  (or  characteristic 
function)  gives  us  a  clear  borderline  between  membership  and  non-membership  for 
classical  set  theory.  In  fuzzy  sets,  the  indicator  function  is  allowed  to  vary  over 
the  range  [0, 1  and  was  retermed  by  Zadeh  as  ’"the  membership  function”  23  In 
this  wa\,  the  uncertainty  as  to  whether  an  object  belongs  to  a  given  set  or  not  is 
expressed. 

If  the  membership  level  is  1,  the  element  or  object  is  definitely  a  member  of 
the  set.  If  the  membership  is  0,  then  it  is  definitely  not  a  member.  However,  if  the 
membership  is  an  intermediate  value,  between  0  and  1.  it  indicates  the  degree  of 
confidence  that  the  object  is  a  member  of  the  set.  A  fuzzy  set,  therefore,  is  a  class 
that  admits  the  possibility  of  partial  membership.  Let  z  —  {z}  denote  a  space  of 
objects.  Then,  a  fuzzy  set  A  in  A'  is  a  set  of  ordered  pairs 

A  =  {z,/ia(x)}  ztX 

where  Va{x)  is  termed  ’’the  grade  of  membership  z  to  A”. 
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Figure  3:  Locality  Index  Membership  Functions 


5.2  Application  Example 

The  fuzzy  theory  will  be  applied  to  the  design  reliability  problem  using  the  cost 
function  developed  using  regression  analysis  on  the  building  variables.  Five  vari¬ 
ables.  including  locality  index  (L).  price  index  (P),  building  type  index  (T).  quality 
index  (Q)  and  technology  index  (TE)  will  be  considered  as  fuzzy  variables,  with 
the  membership  functions  shown  in  Figs.  3-7.  These  membership  functions  are 
assumed  to  be  given. 


5.3  Extension  Principle 

To  find  the  membership  function  the  extension  principle  will  be  used,  as  defined 
below  [lOj. 

Let  X  be  a  Cartesian  product  of  universes,  A'  =  zj  x  ....  x  zr,  and  A\ , . ,  .4, 

be  r  fuzzy  sets  in  Xi,....,Xr  respectively.  The  Cartesian  product  of  .4i,...,.4,  is 
defined  as 

•4i  x .  x  Ar  =  f  r))i(z,, ....,  zr))  (19) 

J  X\  X  ...  XX, 

Let  /  be  a  mapping  from  A'i  >  ....  >  A',  to  universe  Y  such  that  y  -  f{i\ ,  *,). 
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Figure  4:  Price  Index  Membership  Functions 
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Figure  5:  Type  Index  Membership  Functions 
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The  extension  principle  allows  us  to  induce  from  r  fuzzy  sets  .4,  a  fuzzy  set  B  on 
V  through  /.  such  that 

BB{y)  -  sup  [min{fxAl{x -J . ^t(r,))'  (20) 

HB{y)  =  0  if  f~l  =  o 

where  f~l(y)  is  the  inverse  image  of  y  and  o  denotes  an  empty  set.  /Jfi(y)  is 

the  greatest  among  the  membership  values  \  .  >.  fiAr( A'j. . A'r)  of  the 

realization  of  y  using  r-tuples  (x^ . xr). 

As  an  example,  consider  a  variable  combination  1  (see  Table  8)  with  L  = 
,9,P  =  2.76.  T  =  2.95.  H  =  40.  Q  =  2  and  TE  -  1.0.  and  consider  variable 
membership  functions  using  figures  3-7.  (e  g.,  good  locality  index,  low  price  index, 
low  type  index,  low  quality  index  and  low  technology  index),  the  membership 
values  are  m  =  .85, /if  =  .15,  /17  =  ,22./i//  =  l.  (non-fuzzy  variable),  pg  =  .75. 
and  /iff  =  1.  Based  on  the  extension  principle.  Me;  .9.2.76.2.95.40.2.1.01  =  min  ( -So .  .15. 
.22.  1..  .75. 1.)  =  .15  .  The  corresponding  cost  is  calculated  using  the  cost  function 
Eq.  (18).  That  is.  C  =  /(I.  P.T.H.Q.TE )  =  /(.9. 2.76.2.95. 40. 2. 1.0)  =  37.7. 

Each  design  variahle  has  three  possible  membership  functions.  Cost  function  is 
a  function  of  these  design  variables.  Therefore,  different  combinations  of  the  design 
membership  functions  will  result  in  different  cost  membership  functions.  In  this 
example,  three  design  membership  combinations  were  considered:  (a)  good  locality 
index,  low  price  index,  simple  type  index,  low  quality  index,  and  low  technology 
index,  (b)  average  indices  for  all  five  variables,  (c)  Door  locality  index,  high  price 
index,  complicated  type  index,  high  quality  index,  and  high  technology  index 
These  three  combinations  result  in  three  different  cost  membership  functions  as 
shown  in  Table  9.  The  membership  functions  me  (se^  Fig  8),  hmc  (see  Fig.  9) 
HHC  (se?  Fig.  10),  correspond  to  combinations  (a),  (b)  and  (c),  respectively. 

6  Fuzzy  Probability  Approach 

The  cost  function  is  used  to  find  the  estimated  cost  for  each  set  of  design  variables 
.  The  data  for  these  calculations  are  shown  in  Table  9.  with  the  cost  membership 
value  of  all  three  combinations.  Based  on  these  data,  a  probability  mass  function 
is  found  and  shown  in  Figure  11.  From  the  mass  probability  function  and  the 
event  membership  function,  the  probability  of  a  fuzzy  event  can  be  found  as: 


p(A)  =  Y.  V* P » 


(21) 


Figure  10:  High  Cost  Membership  Function 
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Figure  11:  Cost  Probability  Mass  Function 


Accordingly. 

the  probability  of  low  cost  event  is 


P{CL)  =  .0263(  .62  -  .7  -r  .45  -r  .5  -  .5  -  .17  4-  .5  -  .3  -  .22  -  .22  +  .3  -r 
.17  -  .17  -  .15  -  .25  -  .17  -  .04  -  .22  -  .04  -  .15  -  .04)  -  .05263(  62  -  .6  - 
.22  -  .17  -  .22  -  .22)  -  .0789  x  .6  -  .105  x  .17  =  0.3145  (22) 

the  probabilit  y  of  medium  cost  event  is 


P(CM)  =  .0263( .15  -  .6  -  .15  -  .6  -  .45  -t-  .45  -  .45  -  .62  -  .35  - 

.32  -  .3  -  .32  -  .32  -  .15  -  .32  -  .15  -  .45  J-  .15)  -  05263( .15  - 

.6  -  .32  -  .45)  -  .0789  x  .62  -  .105  x  .35  =  .3681  (23) 

and  the  probability  of  high  cost  event  is 

P{CH)  =  0263(  07  -  .07  -  .32  -  .3  -  .32  -  .32  -  35  -  .23  -  .32 

-.4  -  .35  -  .35  -  .45  -  .35  -  .3  ~  .35  -  .32  -  .5)  -  .05263 

(.07  -  .32  -  .35  -  .35)  -  .0789  x  .35  -  .105  x  .35  =  .271  (24) 

The  probability  of  cost  C,  can  be  found  by 

P{C,)=^~  (25) 

L.  Ml 

For  example,  if  low  cost  event  is  considered,  probability  of  C  =  40  will  be: 

P(C  .  40)  =  ^  =  3!  =  0259  (26) 

If  all  cost  ranges  are  considered.  Figure  12  shows  fuzzy  probabilities  for  cost 
values  for  the  case  where  all  variables  have  average  membership  values. 

7  Sensitivity  Analysis 

Sensitivity  analysis  is  used  to  calculate  the  sensitivity  of  the  cost  membership 
function  and  cost  fuzzy  probability  to  change  in  the  design  membership  functions. 
The  cost  is  a  function  of  several  variables.  I'.  Cost  may  be  sensitive  to  change 
in  one  variable  but  not  to  change  in  another.  In  this  section,  a  studv  of  the 


sensitivity  of  the  cost  membership  function  and  fuzzy  probability  to  change  in 
the  independent  design  membership  functions  is  presented  Ten  different  cases  of 
design  membership  functions  are  considered.  These  cases  are  as  follows. 

1.  Good  locality  index  (L) 

2.  Low  price  index  (P) 

3.  Simple  type  index  (T) 

4.  Low  quabty  index  (0) 

5.  Low  technology  index  (TE) 

6.  Poor  locality  index  (L) 

7.  High  price  index  (P) 

8.  Complicated  type  index  (T) 

9.  High  quality  index  (Q) 

10.  High  technology  index  (TE) 

The  ten  cases  mentioned  above  are  '■hosen  to  note  the  effect  that  independent 
design  variables  will  have  on  the  cost  membership  functions  and  on  the  fuzzy 
probability;  for  example  as  they  change  from  average  to  low  or  high  values  (case 
2.  above)  or  (case  7.  above),  respectively. 

Figures  9  and  12  represent  cost  membership  function,  and  fuzzy  probability, 
respectively,  when  all  variables  are  at  their  average  values.  Appendix  B  contains 
the  Output  Data  for  fuzzy  probability  and  membership  values  for  the  average  case 

To  investigate  the  effects  that  the  ues.b»i  rm.*nbe»ship  fu^uoi;  change  will 
have  on  the  cost  membership  function,  a  comparison  is  presented  in  the  following 
section. 


7.1  Analysis  of  the  Results 

Figures  13  and  23  show  the  cost  fuzzy  probability  and  the  cost  membership  func¬ 
tion,  respectively,  that  correspond  to  case  1 .  good  locality  index  These  two  Figures 
show  a  slight  decrease  in  one  intermediate  value,  but  no  increase  in  the  low  cost 
or  decrease  in  the  high  cost  membership  values.  Figures  18  and  28  are  the  cost 
fuzzy  probability  and  the  cost  membership  function  respectively,  that  correspond 
to  poor  locality  index  (cast  6,  above).  This  case  has  no  effect  at  all  on  the  final 
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cost.  Figures  1-1  and  24  are  the  cost  fuzzy  probability  and  cost  membership  func¬ 
tion.  respectively  that  correspond  to  low  price  index(case  2.  above)  Price  index, 
as  shown  in  the  Figures. has  more  effect  than  any  other  variable.  The  study  of 
Figures  14  and  24  and  the  Output  Data  (Appendix  C).  corresponding  to  low  price 
index  show,  for  example,  that  the  membership  value  of  a  low  cost  (e  g..  cost  =  $21) 
was  increased  from  0.32  to  0.67.  and  fuzzv  probability  increased  from  0.042  to 
0  09544 :  also,  the  membership  value  of  a  high  cost(e  g  .  cost  =568)  was  decreased 
from  0.32  to  0.15.  and  fuzzy  probability  was  decreased  from  0.0394  tc  0.02136  On 
the  other  hand  Figures  19  and  29  are  the  cost  fuzzy  probabihtv  and  membership 
function,  respectively,  that  correspond  to  high  price  index  (case  7.  above).  This 
case  has  a  high  effect  on  the  cost  membership  function.  Figure  29  and  the  Output 
Data  (Appendix  C)  show  the  effect  of  high  price  index  on  the  final  cost  mem¬ 
bership  function.  For  example,  the  membership  value  of  low  costfe.g..  cost  =  S2 1 ) 
was  decreased  from  0.32  at  the  average  to  015  at  the  high  price  index,  and  fuzzy 
probability  was  decreased  from  0.042  to  0.01976. 

These  two  examples  show  that  the  effects  of  the  design  variables  are  not  the 
same,  consequentlv.  the  treatment  of  these  variables  should  not  be  the  same  The 
first  variable,  locality  index,  has  less  effect  than  the  the  second  one.  price  index, 
which  requires  more  attention  than  locality  index.  Another  important  point  worth 
mentioning  here  is  that,  what  is  true  for  this  particular  budding  example  might 
not  be  true  for  another  project,  meaning  that  if  locality  index  is  not  critical  m  this 
project  it  may  be  critical  in  another  one  Another  example  of  different  effect  is  a 
variable,  which  has  almost  the  same  effect  on  all  variable  combination  values  .  for 
example  technology  index.  This  variable,  as  shown  in  Figure  22;high  technology 
index'  wih  decrease  almost  all  cost  fuzzy  probability  values  except  tne  high  values, 
which  were  slightly  increased  (e  g.,  cost  =  $68 )  wh«.,e  the  fuzzy  probability  was 
decreased  from  0.C39  to  0  0467  Based  on  sensitive  v  analvsis  the  designer  will 
identify  the  most  critical  variable.  The  designer  will  give  both  emphasis  and 
attention  to  this  based  on  his  interest  and  the  amount  of  risk  he  is  willing  to 
take  For  example,  if  he  is  conservative,  he  will  emphasize  those  variables  that 
decrease  high-cost  membership  and  fuzzy  probability  values,  and  vice  versa 
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Figure  14:  Fuzzy  Probability  of  Cost  Corresponding  to  Lo*  Price  Index(P) 
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Figure  15:  Fuzzy  Probability  of  Cost  Corresponding  to  Simple  Type  Index(T ) 


Figure  18:  Fuzzy  Probability  of  Cost  Corresponding  to  Poor  Locality  In- 
dex(L) 
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Figure  19:  Fuzzy  Probability  of  Cost  Corresponding  to  High  Price  Index(P) 


Figure  20:  Fuzzy  Probability  of  Cost  Corresponding  to  Complicated  Type 
Index(T) 
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Table  9:  Cost  Membership  Functions  and  Probabilities 


Cost  111 
dollars 

square 

foot 

Number  of 
observations 

No. 

Prob  ability 

P 

Membership 

Functions 

LYM1  :  MV M2 

ULC  !  HMC 

HVM 3 

VHC 

Variable 

Combinations 

10 

1 

0.02631 

0.62 

0.00 

0.00 

37 

11 

1 

0.02631 

0.70 

0  00 

0.00 

35 

14 

1 

0.02631 

0.45 

0.15 

0.07 

38 

15 

2 

0.05263 

0.62 

0.15  i 

0.07 

36.27 

17 

1 

0.02631 

0.10 

0.00  ; 

0.00 

28 

20 

1 

0.02631 

0.50 

0.00  ; 

0.00 

25 

21 

2 

0.05263 

0.60 

0.32  i 

0.15 

24.31 

23 

3 

0 07890 

0.60 

0.62  ■ 

0.32 

6.10.32 

24 

2 

0.05263 

0.22 

0.00 

0.00 

2.22 

27 

1 

0.02631 

0.17 

0.00  : 

0  00 

34 

31 

1 

0.02631 

0.50 

0.60 

0  30 

33 

32 

1 

0.02631 

0.22 

0.45 

0.32 

4 

33 

1 

0.02631 

0  22 

0.45 

0.32 

18 

34 

1 

0.02631 

0.22 

0.45  , 

0.35 

26 

36 

4 

0.10526 

0  17 

0.35 

0.35 

16.17.21.30 

37 

1 

0.02631 

0.30 

0.62 

0.35 

5 

38 

2 

0.05263 

0.17 

0.32 

0.35 

1.9 

39 

1 

0.02631 

0.17 

0.35 

0.32 

29 

40 

2 

0.05263 

0.22 

0.45  ' 

0.35 

7.19 

41 

1 

0  02631 

0.17 

0.32 

0  40 

23 

43 

1 

0.02631 

0.15 

0.30 

0.35 

11 

45 

1 

0  02631 

0.25 

0.32  ! 

0.35 

12 

46 

1 

0.02631 

0.17 

0.32  | 

0.45 

15 

48 

1 

0.02631 

0.04 

0.15  ; 

0.35 

20 

56 

1 

0.02631 

0.22 

0.32 

0.30 

8 

62 

1 

0.02631 

0.04 

0.15 

0.35 

13 

68 

1 

0.02631 

0.15 

0.30  : 

0.53 

3 

86 

1 

0.02631 

0.04 

0.15 

0  50 

14 

|  1-Low  Variable  Membership  Combination 
|  2-Medium  Variable  Membership  Combination 
3-High  Variable  Membership  Combination 


8  Conclusion 


The  actual  cost  deviation  from  the  feasibility  estimate  is  largely  due  to  a  lack  f 
information  concerning  the  project  and  its  many  design  variables.  This  problem 
demonstrates  the  necessity  of  having  a  model  that  can  handle  it  and  devise  a  model 
that  can  derive  more  accurate  results  at  this  stage.  Regression  analysis  is  used  to 
find  the  relation  between  historical  data  and  design  variables.  This  relation  can  be 
helpful  if  we  know  the  values  of  the  independent  variables  precisely,  however,  in  the 
real  world  this  is  not  always  the  case.  In  the  process  of  design,  there  are  many  cases 
in  which  the  available  information  concerning  the  value  of  the  variable  in  the  future 
will  be  known  only  by  some  verbal  expression,  such  as  ’’high,”  ’Tow,’'  ’’good,” 
or  ’’bad.”  These  expressions  cannot  be  used  directly  in  such  a  function.  Fuzzy 
set  theory  is  applied  by  assigning  membership  functions  to  these  design  variables 
according  to  design  experience  and  knowledge  concerning  them.  The  cost  will  be 
found  using  the  cost  function,  and  using  extension  principle,  the  cost  membership 
value  will  be  assigned  based  on  the  independent  variable  membership  function. 
The  example  given  was  that  of  a  building,  but  this  methodology  can  be  applied 
to  any  kind  of  project  as  long  as  the  design  variables  are  known  and  the  data 
for  dependent  and  independent  variables  are  available.  A  computer  program  was 
developed  to  find  the  cost  function,  cost  probability  mass  function  and  membership 
function  (see  Appendix  A).  The  cost  function  for  any  project  can  be  found  with 
both  the  same  program  and  procedure. 


8.1  Required  Information  to  Use  the  Design  Relia¬ 
bility  Model 

The  developed  method  to  estimate  the  project  cost  at  the  feasibility  stage  is  not 
limited  to  a  building  projects,  as  given  in  the  report,  but  can  be  applied  to  any 
type  of  project,  including  pile-driving  projects  if  the  following  data  are  provided: 

1.  Dependent  variable  (cost  variable),  this  variable  can  represent  cost  per  pile, 
cost  per  square  foot,  etc. 

2.  Independent  design  variables  that  affect  the  cost  (\\,i  =  1. . m):  the  value 

of  these  variables  depend  on  cost  variable  definition  and  how  they  will  con¬ 
tribute  in  that  cost  specified  in  (1). 

3.  Historical  data  for  such  design  variables. 

4.  Future  prediction  of  variables. 
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function 


Figure  29:  Membership  Function  of  Cost  Corresponding  to  High  Price  In 
dex(P) 
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Figure  30:  Membership  Function  of  Cost  Corresponding  to  Complicated 
Type  Index(T) 


The  last  item  is  needed  to  develop  membership  functions  of  the  variables.  It 
is  enough  to  know  verbal  information,  such  as  "good."  ’’bad.”  ’’low,”  and  so  on. 
It  is  also  helpful  to  give  a  range  of  variables  for  maximum  and  minimum,  to  help 
assign  membership  values  for  such  variables  within  that  range. 

To  verify  the  model,  a  historical  of  actual  and  predicted  design  variables  of 
past  projects  will  be  of  great  help. 
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APPENDIX  A 
Program  Listing 

................  hxis  PROGRAM  *•••*«*••••••«*»••*•••••••**•••• 

DIMENSION  A (50, 50) ,C (50) ,X(50) ,  XP (50, 50) ,  YP  (50) , SUM (50) , SUM1  (50) 
DIMENSION  SUM2 (50) , SUM3 (50) , EMU (50, 50) , CMU(50> , CP (50) , yl (50) 
INTEGER  0(50) ,CO(50) 

INTEGER  H,  Uni,  i 2,  i3 


character*40  title 
character*12  :nfn 
character*12  outfnl 
character*12  outfn2 
character*12  out£n3 

C  character*12  outfn4 

••*.««The  program  will  prompt  user  to  input  the  following  files 

*  1-Input  data  file 

*  2-file  that  contains  regression  analysis  output  which  includes 

*  3-file  that  contains  cost  membership  function  and  cost  PMF 

*  and  graphs  of  the  data 

*  4-file  that  contains  input  data 


print*, 'please  enter  the  output  title  max  40  characters' 
read(«, 905) title 
905  format (A40) 

700  print*, 'please  enter  input  file' 
read (*, 900 , ERR-700) infn 
900  format (A12) 

OPEN (UNIT-50, riLE-infn, STATUS-' OLD' ) 

print*, ' please  enter  regression  analysis  output  file  ' 
read(*, 900) outfnl 

print*, 'please  enter  probabilities  and  membership  functions 
.output  file  ' 
read (* , 900) outfn2 
C 

print", ' please  enter  data  output  file  ' 
read(*, 900) outfn3 

C  print*, ' please  enter  data  output  file  ' 

C  read (*, 900 ) outfn4 

OPEN (UNIT-60, EILE-OUtfnl , STATUS-' UNKNOWN' ) 

OPEN (UNIT-65 , FILE-outf n2 , STATUS-' UNKNOWN' ) 

OPEN (UNIT-66, riLE-outfn3, STATUS-' UNKNOWN' ) 

C  OPEN (UNIT-67, FILE-outfn4, STATUS-' UNKNOWN' ) 

C  WRITE v60, 905) title 

WRITE (65, 905) title 
WRITE (66, 905) title 
C  WRITE (67, 305) title 

.......................... INpu  QATA*. 

•Read  number  of  variables  and  number  of  data  points 

C  PRINT*,  '  PLEASE  ENTER  THE  NUMBER  OF  VARIABLES' 

READ  (50,  •)  M 

Cll  PRINT-, ' PLEASE  ENTER  NUMBER  OF  DATA  POINTS' 

11  READ (50, * ) L 

N-M+l 


IF  (L  .  LE  .  M)  THEN 


N  MUST  BE  >M 


PRINT*, ' 

GO  TO  11 
ELSE 

*  *Read  dependent  variable  , independent  variables  and  variables  membership 

*  values 

DO  12  I-1,L 

READ  (50,  *)  TP  (I)  ,  (XP  (J,  I)  ,  J-1,M) 

READ (50, *)  <EMU(J, I) , J-1,M) 

WRITE (66, *)  re  (I) ,  (XP (J, I) , J-1,M) 

WRITE ( 66 , * ) (EMU (J, I) , J*1,M) 

12  CONTINUE 

ENDIF 

•Call  subroutines  to  construct  simultaneus  equation  and  solve  for  coefficients 

CALL  MULTLR (A, XP,YP,M,L,N,C) 

CALL  ORDER (S , O, A, N) 

CALL  ELIM (S, A, 0,N) 

CALL  SOLVE (A, C,X,0,N) 

•Write  the  resulted  cost  function 

C  WRITE (65, 204) 

204  format (2x, 75 ('«') ) 

C  WRITE (65,200) 

C  WRITE (65  252) 

C252  FORMAT (3 OX, 'THE  COST  FUNCTION'/) 

C  WRITE (65, 204) 

C  WRITE (65,200) 

C  WRITE (65, 299) X(l) , (X(J) , (J-l) , J-2,M+1) 

C299  FORMAT ( IX, '  Y-' , F10 . 5, ' +' , 10 ( ' ( ' , F10 . 5, ' ) ' , IX, ' X' , 12 , ' ♦ ' ) ) 

•Call  membership  subroutine  to  find  cost  membereship  function  based  on 
•the  input  independent  variables 

* 

CALL  MEMBERSHIP (XP, CO, X,L,M,N, EMU, CMU) 

•Call  sorting  subroutine  to  sort  the  output  data  to  simplify  grouping 
•subroutine  to  group  the  equal  points  and  have  one  membership  value  for 
•each  point 

CALL  SORTING (CO, CMU,L) 

CALL  GROUPING (CO, CMU,L) 

•Write  the  cost  cofficient  values 


WRITE (60, 200) 

200  format (2x, 75 ('■')/ ) 

WRITE (60, 205) 

205  format (30X, 'COST  FUNCTION  COEFFICIENTS'/) 

WRITE (60, 200) 

WRITE ( 60 , 97) X ( 1 ) 
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97  FORMAT ( 5X, ' a  0=',F10.6/) 

DO  14  I»2,N 

WRITE  (60,  99)  { I  - 1 )  ,  X  (I) 

99  FORMAT (5X,  '  a'  ,  12,  ' ='  ,  F10 . 6/) 

14  CONTINUE 

•Quantify  Error  associated  with  the  selected  linear  function 

ST=0 . 0 
SR=0 . 0 

DO  16  I-l.L 

SUM ( I ) =0 . 0 
SUM1 (I) *0 . 0 

ST-ST+ (YP (I) -YMEAN) **2 
DO  18  J=1 , N 

IF (J.EQ. 1) THEN 
SUM(I)«X(1) 

ELSE 

SUM1  ( I )  “SUM1  (I)  +X  (J)  •  (XP  (  (J-l)  ,  I)  ) 

END  IF 

18  CONTINUE 

SUM2 (I) “SUM ( I ) +  SUM1 (I) 

SR-SR+ (YP (I) -SUM2 (I) ) **2 

16  CONTINUE 

•Write  the  result  of  the  error  quantification 
WRITE (60, 200) 

WRITE  (60, 210 

210  format ( 25x, ' QUANTIFICATION  OF  ERROR  OF  LINEAR  REGRESSION' ! ) 

WRITE (60, 200) 

WRITE (60, 225) ST 

225  format (2x, ' THE  TOTAL  OF  SUM  OF  SQUARES  OF  RESIDUALS  St-',F10.3) 

WRITE (60, 230) SR 

230  FORMAT <2X, ' THE  SUM  ODF  SQUARES  OF  THE  RESIDUALS  Sr*',F10.3) 

SXY-SQRT (SR/ (L- (M+l) ) ) 

R2« (ST-SR) /ST 
R-SQRT (R2) 


WRITE (60, 235) SXY 

235  FORMAT (2X, ' THE  STANDARS  ERROR  OF  ESTIMATE  Sx/y-' , f 10 . 3 i 

WRITE  (60, 240)  R2 

240  FORMAT (2X, ' THE  COEFFICIENT  OF  DETERMINATION  r'2=',F6.4! 


WRITE (60, 245) R 


245  FORMAT (2X,  '  THE  CORRELATION  COEFFICIENT  r-',F6.4) 

WRITE (60, 207) 

207  format  ('  1' 2x,  75  ('-')  ) 

WRITE (60,200) 

WRITE (60, 250) 

250  FORMAT (3 OX, ' THE  FINAL  COST  FUNCTION'/) 

WRITE (60, 204) 

WRITE (60, 200) 

WRITE (60, 98) X(l) , (X(J) , (J-l) , J-2,M+1) 

98  FORMAT (IX, '  Y-'  ,  F10 . 5,  '  +'  ,  10 ('  ('  ,  F1Q . 5,  '  )  '  ,  IX,  '  X' , 12,  ' +'  ) ) 

30  CONTINUE 


STOP 

END 

«•••«*. ••••••••-•••••••MULTILINEAR  REGRESSION  SUBROUTINE 

SUBROUTINE  MULTLR  (A,  XP ,  YP,  M,  L,  Nf  C) 

DIMENSION  A (50, 50) , XP (50, 50) , YP (50) ,C  (50) 

INTEGER  H 

DO  10  I-1,N 
DO  20  J-l, I 

IF  ( 1 . EQ . 1 ) THEN 
A(l, 1) -L 
ELSE 

SUM-0 

DO  30  H— 1 ,  L 
XP  (0 ,  H)  — 1 

SUM— SUM+ (XP (I-1,H) ) * (XP ( J-1,H) ) 

30  CONTINUE 

A ( I , J) —  SUM 
A(J, I) «  SUM 
END  IF 

20  CONTINUE 

SUM-0 . 0 

IF (I .EQ. 1) THEN 

DO  25  H— 1 , L 

SUM-SUM+Y? (H) 

25  CONTINUE 

C (1) -SUM 

ELSE 

DO  40  H-l , L 

SUM-SUM+YP (H) * (XP (I-1,H) ) 

40  CONTINUE 
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C(I)»SUM 


END  IF 
10  CONTINUE 

WRITE (60, 360) 

360  FORMAT (2X, 75 ('«')/) 

WRITE (60, 370) 

370  FORMAT (4 OX, 'OUTPUT  '/) 

WRITE (60, 380) 

380  FORMAT (2X, 75 ('*') /) 

DO  60  I  -1,N 

WRITE (60, 100)  (A(I, J) , J-1,N) ,C  (I) 

100  FORMAT (2X, 10 (2X,F14 . 6) , 2X, '  C-',F14.6) 

60  CONTINUE 

RETURN 

END 


**».««*«*.*0RDER  subroutine 

SUBROUTINE  ORDER (S, O, A, N) 

DIMENSION  A(50, 50) , S (50) 

INTEGER  0(50) 

DO  10  I-1,N 

0(I)-I 

S(I)-ABS  (A  (1,1)) 

DO  20  J«2,N 

IF  (ABS (A (I, J) ) .GT.S (I) )  THEN 
S  (I)  -ABS  (A  ( I ,  J)  ) 

END  IF 

20  CONTINUE 

10  CONTINUE 
RETURN 
END 


****** .......... ••••SUBROUTINE  ELIM**** 

SUBROUTINE  ELIM (S, A, O, N) 

DIMENSION  A (50, 50) ,C(50) ,S(50) ,F(50,50) 
INTEGER  0(50) 

INTEGER* 2  I,J,K 

DO  10  K-1,N-1 

CALL  PIVOT (S,A,0,N,K) 

DO  20  I-K+1,N 

F(0(I)  ,K)-A(0(I)  ,K)  /A  (O  (K)  ,K) 

DO  30  J-K+1,N 
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A(0(I)  , J) -A (0  ( I )  ,J)  -F  (0(1)  ,  K)  *A(0(K)  ,J) 


30  CONTINUE 

20  CONTINUE 

10  CONTINUE 

DO  40  K-1,N-1 

DO  50  I-K+1,N 

A(0(I)  ,K)-F  (0(1)  ,  K) 

50  CONTINUE 

40  CONTINUE 

RETURN 

END 

S0LVE  SUBROUTINE 

SUBROUTINE  SOLVE (A, C, X, O, N) 

DIMENSION  A(50,50) , C (50 ) , X ( 5 0 ) 

INTEGER  0(50) 

X(1)-C(0(1)  ) 

DO  10  1-2, N 

SUM-0 . 0 

DO  20  J-l.I-1 

SUM-SUM+A (0(1) , J) *X(J) 

20  CONTINUE 

X(I)-C(0(I)  )  -SUM 

10  CONTINUE 

X  (N)  —X  (N)  /A  (O  (N)  ,  N) 

DO  30  I-N-l , 1,-1 

SUM-0 . 0 

DO  40  J— 1  +  1 , N 

SUM-SUM+A (0(1) , J) *X ( J) 

40  CONTINUE 

X(t)-(X(I) -SUM)  /A (0(1)  ,  I) 

30  CONTINUE 

RETURN 

END 


SUBROUTINE  PIVOT 
SUBROUTINE  P IVOT  ( S ,  A,  0,  N,  K) 
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DIMENSION  A(50, 50) , S  (50) 

INTEGER  0(50) 

INTEGER* 2  K, PIVIT, II, IDUM 
PIVIT  -  K 

BIG-ABS (A(0(K) ,K) /S (0 (K) ) ) 

DO  10  II«K+1,N 

OUMMY-ABS (A(0(II) ,K) /S (O(II) ) ) 

IF  (DUMMY. GT. BIG)  THEN 

BIG-DUMMY 

PIVIT-II 

END  IF 

10  CONTINUE 

IDUM-0 (PIVIT) 

0 (PIVIT) -O (K) 

0(K)»IDUM 

RETURN 

END 


.......... -SUBROUTINE  MEMBERSHIP  FUNCTION****"************** ******** 

SUBROUTINE  MEMBERSHIP (XP , CO, X, L, M, N, EMU, CMU) 

DIMENSION  XP (50, 50) ,X(50) , SUM3 (50) , COl (50) , DIF (50) , CMU (50) 

DIMENSION  EMU (50, 50) 

INTEGER  CO  (50) 

•Find  the  coat  corresponding  to  each  independent  variable  set  and  find  cost 
•membership  value 

DO  15  1*1,1. 

SUM3 (I) -0 

DO  110  J-2,N 

SUM3 (I) -SUM3 (I) +X (J) «XP ( (J-l)  ,  I) 

110  CONTINUE 

CO (I) -X (1) +SUM3 (I) 

COl (I) -X(l) +SUM3 (I) 

DIF (I) -COl (I) -CO (I) 

CMU ( I ) -EMU ( 1 , I) 

DO  114  J-1,M-1 

IF (CMU (I) .LT.EMU ( ( J+l) , I) ) THEN 
CMU (I) -CMU (I) 

ELSE 

CMU  ( I )  -EMU  (  (J-t-1)  ,  I) 

END  IF 

114  CONTINUE 

IF (DIF ( I ) . GT .0.5) THEN 
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C0(I)-C0(I)+1 


30 


40 


c 

C255 

C 

C257 

C 

C260 

C 

C 

C258 

C 

C256 

C 

C265 

C 

C275 


ELSE 

ENDIF 

CONTINUE 

RETURN 

END 


••••■it************** « SUBROUT INE  SORTING 

SUBROUTINE  SORTING (CO, CMU, L) 

DIMENSION  CO (50)  ,  CMU(50) 

INTEGER  L, M, CO 

LOGICAL  SORTED 

M-L-l 

SORTED-. FALSE. 

IF  ( .NOT. SORTED)  THEN 
SORTED- . TRUE . 

DO  40  1-1 , M 

IF  (CO ( I )  . GT . CO (1+1) ) THEN 
TEMP -CO ( I ) 

TEMPM-CMU (I) 

CO(I)-CO(I+l) 

CMU (I) —CMU (1+1) 

CO (I+l)  -TEMP 
CMU ( I + 1 ) -TEMPM 
SORTED-. FALSE. 

ENDIF 

CONTINUE 

GO  TO  30 


ENDIF 

WRITE (65, 255) 

FORMAT  <'1'2X,  75  ('-')) 

WRITE (65, 257) 

FORMAT (2X,  75 ('-')/) 

WRITE (65, 260) 

FORMAT (I5X, 'THE  COSTS  AND  THEIR  CORRESPONDING  MEMBERSHIP 
1  FUNCTIONS'/) 

WRITE (65,258) 

FORMAT (2X, 75 ('-') ) 

WRITE (65, 256) 

FORMAT (2X, 75 ('-')//) 

WRITE (63, 263) 

FORMAT ( 12X, '  COST' , 2X, '  | ',1X, ' MEMBERSHIP  VALUES  ! ' /) 

WRITE (65, 275) 

FORMAT  ( 2'A,  75  ('  -'  )  ) 
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c 


DO  60  1-1,1. 


C  WRITE (65, 270) CO (I) , CMU (I) 

C270  FORMAT (9X, 1 6 , 4X, '  I ' , 6X, F7 . 5 , 1  OX,  '  I  '  ) 

C  WRITE (65, 275) 

C60  CONTINUE 

RETURN 

END 


**•»«**»«***»»«»*»********•* SUBROUTINE  GROUPING***** 

SUBROUTINE  GROUPING (CO, CMU, L) 

DIMENSION  CO (50) , CMU (50) ,NUM(50) , AP (50)  ,  CP (50) 
INTEGER  CO, L, K, K1 , II , NUM, line 
CHARACTER* 60  GTIT 

K-L-l 

Kl-L-1 

11-1 

DO  20  I-1,K 

CO(I)-CO(Il) 

NUM(I)  -1 

AP (I) -REAL (NUM (I) ) /REAL (L) 

CMU (I) -CMU (II) 

11  IF (CO  (II)  . EQ . CO (11+1) ) THEN 

CO(I)-COdl) 

NUM ( I ) —NUM ( I ) +1 

AP (I) -REAL (NUM (I) ) /REAL (L) 

IF  (CMU  (ID  .GT. CMU  (11  +  1)  )  THEN 

CMU  (I) -CMU  (ID 


ELSE 

CMU (I) -CMU (11+1) 
END  IF 
Il-Il+l 

C  K-K-l 

Kl-Kl-1 

IF (II .EQ. L) THEN 

GO  TO  21 

ELSE 
END  IF 
GO  TO  11 

ELSE 

IF  ( 1 1 . EQ . L) THEN 
GO  TO  21 


ELSE 

ENDIF 


ENDIF 


55 


IF (Il.EQ.L)THEN 
GO  TO  21 
ELSE 
END  IF 
Il-Il+l 

20  CONTINUE 

21  CONTINUE 

C  WRITE {67, *)' cost  membeship  values' 

C  do  407  i-1 , kl  +  1 

C  WRITE (67, *) co (i) , emu (i) 

C407  continue 

c****This  part  will  find  probability  of  Ci's 
c 

C  CP-COST  PROBABILITY 

C  TCMU  -TOTAL  SUM  OF  MEMBERSHIP  FUNCTIONS  FOR  COST  EVENT 
TCMU-0 . 0 
DO  400  1-1 , Kl+1 

a  . 

TCMU -TCMU +CMU*(  I ) 

400  CONTINUE 

C  WRITE (65, *)' total-' , tcrau 

C  WRITE (67, *)' cost  and  coat  furry  probability' 

DO  402  I-1,K1+1 

CP  (I)-CMU(I) /TCMU 

C  WRITE ( 67 , * ) CO ( I ) , CP ( I ) 

402  CONTINUE 

•  Prompt  the  user  to  enter  line  length  , this  is  required  for  the  plot  of  the 
*data  this  length  will  be  the  maximumum  length  in  the  graph  and  all  the  data 

*  will  be  scaled  according  to  that  length. 


C  write (67, *)' cost  probability  ' 

•call  subroutine  graph  to  plot  the  cost  vs  their  probabilities 
WRITE (65, 280) 

200  FORMAT (2X, 75 ('-') ) 

WRITE (65, 290) 

290  FORMAT (3 OX, ' THE  FINAL  RESULTS') 

WRITE (65, 289) 

289  FORMAT (2X, 75 ('-')) 

WRITE (65, 293) 

293  FORMAT (2X, 75 ('-') ) 

WRITE (65, 295) 

295  FORMAT(5X, 'COST' ,2X, ' | ' , IX, 'No.  OF  OBS . ' , IX, ' I ' , 3X, ' PMF' , 3X, ' | ' 
1, IX, 'MEMBERSHIP  VALUES' ,' I ’,' FUZZY  PROBABILITIES') 

WRITE (65, 293) 

DO  30  1-1, Kl+1 

WRITE (65, 298) CO (I) , NUM (I) , AP (I) , CMU(I)  ,CP (I) 

298  FORMAT (5X, 14, IX,' |' ,5X,I3,5X, ' ! ',1X, F7 . 6, IX, • | ' , 5X,F7 . 6, 6X, 

1' I ' , F7 . 6) 
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o  o 


WRITS (65, 293) 


30  CONTINUE 

C  WRITE (67, •)' PMF  DATA  ' 

do  409  i-l.kli-l 

WRITE  (67,  *)  co  (i)  ,  ap  d) 
C409  continue 


print*,'  Plot  Option  ' 

print* ,' Please  enter  (1)  if  plot  is  needed' 

print*,'  or (2)  if  plot  is  not  needed' 

read* , i2 
if { i 2 . eq , 1 ) then 

print*, 'Line  max  length  Option' 

print *,' Please  enter  (1)  To  use  default  (Max.  L«40) ' 

print*,'  (2)  To  use  new  value  for  Max.  L' 

read* , i3 

if (i3 . eq. 1) then 
line*40 
else 

if  ( i 3 . eq. 2) then 

print* ,' Please  enter  line  max.  length' 
read" , line 
else 

go  to  411 
endif 
endif 
GTIT-' PMF' 

call  graph  (co,  ap,  kl  +  1,  line,  gtit,  i3) 


GTIT-' Cost  Membership  Function' 
call  graph  (co,  emu,  kl  +  1,  line,  gtit,  i3) 

GTIT»' Cost  Fuzzy  Probabilities' 
call  graph (co, cp, kl  +  1 , line, GTIT, i3) 
else 

if ( i2 . eq . 2 ) then 
RETURN 

else 

go  to  414 
endif 
endif 

return 

END 


*  This  subroutine  prints  a  bar  graph  using  an  array  of 
*N  elements  with  maximum  line  size 

subroutine  graph (co, value, n, line, gtit, i 3' 

integer  n,  line , i , k , co ( 50 )  .nicest,  i3 
real  value ( n) , max , scale , y 1 ( 3 0 i 
logical  r.egnum 
character *20C  bar 
CHARACTER* 60  GTIT 

*D***A****R****E****E****N** 

* 

*  Find  maximum  and  check  for  error  condition. 


•p^ot 

414 

411 
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if  ( line . gt . 200 . or . line . It . 1) then 
print*, 'Line  length  error', line 
else 

negnum- .false . 
max-value  ( 1 ) 
do  10  i-l,n 

if(value(i) .lt.0.0)then 
negnua* . true . 

else 

If (value (i) . gt .max) max»value (i) 
end  if 

13  continue 

* 

*  Fill  bar  with  dashs 

* 

if ( . not . negnum) then 
do  20  i=l , line 
bar  <i :  i) 

20  continue 

* 

*  Scale  data  values  and  print  bar. 

* 

scale-real (line) /max 
WRITE ( 65 , 22 ) GTIT 

22  format (' 1' 30x, 'plot  of ' , lx, A60//) 

write (65 , 23) gtit 

23  format ( 3Cx, A60 ) 
if { i 3 . eq . 1 ) then 

yl ( 1 ) -real (max) /5 
yl  (2) -yl  (1) *2 
yl (3) -yl (1) *3 
yl ( 4 ) -yl (1) *4 
yl  (S)  --max 

WRITE ,65,24) (yl ( i) ,i-i,5) 

24  format (20x, 5 (4x, f4 . 3) ) 

WRITE (65, 26) 

2  6  format  llSx,  '  I - ! - I - I - 

else 

endif 

*  to  print  xtitle 

mccst-n/2 
do  30  i-l, n 

k*nmt  (value  (i)  'scale) 
if (i . eq.mcost) then 

WRITE (65, *)' cost ', co  (i) , '  1', bar (1 : k) 

else 

WRITE  ( 65 ,  * )  '  '  ,  co  ( i )  ,  '  |  '  ,  bar  ( 1 :  k) 

endif 
continue 
endif 
endif 
return 
end 
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APPENDIX  B 

Output  Data  Corresponding  to  Low.  Average,  and  High  Design 
Variable  Membership  Combinations 
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OUTPUT 


38 . 000000 

40 . 919991 

87.049995 

100.539986 

230 . oc::: 

0 

94 . 000000 

41 . 350002 

1288 .060059 

40 . 919991 

44 . 281006 

93 . 509293 

108.509201 

294 . 3300- 

€ 

101 .269989 

44 . 457493 

1388 . 371704 

S'7  .  049995 

93 . 509293 

205 .269714 

231 . 988861 

693 . 93011 

219.219986 

95 . 615997 

3078  .  763672 

100.539986 

108.509201 

231 . 988861 

278.519623 

806 . 48  999 

0 

259.260040 

110 . 477974 

3611.197266 

280 . 000000 

294 . 390076 

691  .  930115 

806.489990 

4746. ::::: 

0 

697 . 000000 

297 . 149994 

10626.958984 

94 . 000000 

101.269989 

219 . 219986 

259.260040 

697.00000 

0 

280 . 000000 

107 . 500000 

3699 . 840332 

41 . 350002 

44 . 457493 

95 . 615997 

110.477974 

297.14999 

4 

10" .500000 

47.762501 

1527 . 611084 

COST 

FUNCTION  COEFFICIENTS 

tVKsv  =  xsfli  =  =  =  m  =  = 

======= 

a  0—84 . 976936 
a  1-  26 . 670523 
a  2-  11.234761 
a  3-  6.235301 
a  4-  0.201219 
a  5-  5.039278 
a  6-  31.219091 


QUANTIFICATION  OF  ERROR  OF  LINEAR  REGRESSION 


THE 

THE 

THE 

THE 

THE 


TOTAL  OF  SUM  OF  SQUARES  OF  RESIDUALS  St- 


SUM  ODF  SQUARES  OF  THE  RESIDUALS  Sr- 

STANDARS  ERROR  OF  ESTIMATE  Sx/y- 

COEFFICIENT  OF  DETERMINATION  r*2«0 

CORRELATION  COEFFICIENT  r-0 


53183.535 
30.700 
0  . 

.  9994 
.  9997 


995 


MJ 


Y*  -84 . 97694+  (  26.67052)  X  1+ (  11.23476)  X  2+ (  6.23530)  X  3+ ( 

4  + {  5.03928)  X  5+ (  31.21909)  X  6  ' 


0 . 20122 
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OUTPUT  CORRESPONDING  TO  VARIABLES  AT,  LOW  yfXKBILVLSKTP^VTiZ.afA 


COST 

I  No.  or  OBS 

.  1  PMF  | 

MEMBERSHIP 

VALUES  1  FUZZY  PROBABILITIES 

10 

1 

1  .026316  | 

.620000 

1 .077500 

11 

1 

1  .026316  | 

.700000 

1 .087500 

14 

— 

1 

I  .026316  | 

.450000 

1.056250 

IS 

2 

!  .052632  1 

.620000 

1,077500 

17 

1 

1  .026316  I 

.100000 

1 .012500 

20 

1 

I  .026316  1 

.500000 

1 .062500 

21 

. 2 

1  . 05J2632  I' 

.600000 

1 .075000 

23 

3 

1  .078947  | 

. 600000 

1.075000 

24 

2 

I  .052632  | 

.220000 

1.027500 

27 

1 

I  .026316  I 

.170000 

I  ,021250 

31 

1 

I  .026316  | 

.500000 

1.062500 

32 

1 

I  .026316  1 

.220000 

1.027500 

33 

1 

1  .026316  I 

.220000 

1  .027500 

34 

1 

1  .026316  I 

.220000 

I  .027500 

36 

4 

I  .105263  I 

.170000 

1.021250 

37 

1 

I  .026316  | 

.300000 

1.037500 

38 

2 

1  .052632  I 

.170000 

I .021250 

39 

1 

|  .026316  | 

.170000 

1 .021250 

40 

2 

1  .052632  I 

.220000 

1.027500 

41 

1 

1  .026316  | 

.170000 

1.021250 

43 

1 

1  .026316  I 

.150000 

1 .018750 

45 

1 

I  .026316  | 

.250000 

1 .031250 

46 

1 

I  .026316  1 

.170000 

1.021250 

48 

1 

I  .026316  1 

.040000 

I .005000 

56 

1 

1  .026316  | 

.220000 

1 .027500 

62 

1 

1  .026316  I 

.040000 

I .005000 

68 

1 

I  .026316  I 

.150000 

I  .018750 

86 

1 

1  .026316  | 

.040000 

1.005000 

THE  FINAL  RESULTS 
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cost 


PMF 

.021  .042  .063  .084  .105 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


63 


Membership  Function 


Fuzzy  Probabilities 


coat 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


.018  .035  .053  .070  .088 
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PU7PU7  CORRESPONDING  TO  VARIABLES  AT  THEIR  MEDIUM  MEMBERSHIP  VALUES 


THE  FINAL  RESULTS 


COST 

1  No. 

OF  OBS 

.  1  PMF  1 

MEMBERSHIP 

VALUES  1  FUZZY  PROBABILITIES 

10 

1 

1 

1  .026316  1 

. 000000 

1  .  000000 

11 

1 

1 

1  .026316  1 

.000000 

1  . 000000 

14 

1  — 

1 

1  .026316  1 

.150000 

1  .  019711 

15 

1 

2 

1  .052632  | 

.150000 

1  .019711 

17 

1 

1 

1  .026316  1 

.000000 

1  .000000 

20 

1 

1 

1  .026316  1 

.000000 

1  .  000000 

21 

1 

2 

1  .052632  I 

. 32CC0C 

i  .  042050 

23 

1 

3 

1  .078947  | 

. 620000 

I  . 081472 

24 

1 

2 

1  .052632  | 

.000000 

1  .  000000 

27 

1 

1 

I  .026316  ! 

.000000 

1 . 000000 

31 

1 

1 

I  .026316  I 

.600000 

1  . 076644 

32 

1 

1 

1  .026316  I 

.450000 

1  .059133 

33 

1 

1 

1  .026316  1 

.450000 

1  .  059133 

34 

1 

1 

1  .026316  | 

.450000 

1  .  059133 

36 

1 

4 

1  .105263  1 

.350000 

I  .045992 

37 

1 

1 

1  .026316  1 

.620000 

1  .081472 

38 

1 

2 

1  .052632  1 

.320000 

i  .042050 

39 

1 

1 

1  .026316  I 

.350000 

1  .  C45992 

4C 

1 

2 

1  .052632  l 

. 450000 

1  .  059133 

41 

1 

1 

I  .026316  I 

.320000 

1  .  042050 

43 

1 

1 

I  .026316  1 

. 300000 

1  .  039422 

45 

1 

1 

1  .026316  1 

. 320000 

1  . 042050 

46 

1 

1 

1  .026316  ; 

. 320000 

1  .  042050 

48 

1 

1 

1  .026316  I 

. 1500C0 

1  .  019711 

56 

1 

1 

1  .026316  1 

. 320000 

1  . 042050 

62 

1 

1 

1  .026316  1 

. 150000 

i  .019711 

68 

1 

1 

1  .026316  I 

. 300000 

1  .  039422 

86 

1 

1 

1  .026316  1 

.  150000 

1  .019711 

66 


PMF 


cost 


.021  .042  .063  .084  .105 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 
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Membership  Function 


cost 


•124  .248  .372  .496  .620 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


68 


Fuzzy  Probabilities 


cost 


.016  .033  .049  .065  .081 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 
45 

48 

56 

62 

68 

86 
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OUTPUT  CORRESPONDING  TO  VARIABLES  AT  THEIR  HIGH  MEMBERSHIP  VALUES 


THE  FINAL  RESULTS 


COST 

1  No. 

OF  OBS 

.  |  PMF  I 

MEMBERSHIP 

VALUES  I  FUZZY  PROBABILITIES 

10 

1 

1 

1  .026316  | 

. 000000 

1 . ooocoo 

11 

1 

1 

I  .026316  I 

. 000000 

1 . 000000 

14 

1  — 

1 

1  .026316  i 

. 070000 

1  .010219 

15 

1 

2 

1  .052632  | 

. 070000 

1  .  010219 

17 

1 

1 

1  .026316  | 

.  000000 

1  .000000 

20 

1 

1 

1  .026316  I 

. 000000 

1  .  000000 

21 

1 

2 

1  .052632  1 

.150000 

1  .  021898 

23 

1 

3 

1  .078947  | 

320000 

1  .046715 

24 

1 

2 

1  .052632  | 

. 000000 

I  . oooooo 

27 

1 

1 

1  .026316  1 

.000000 

1  .  oooooo 

31 

1 

1 

1  .026316  I 

.300000 

1  .043796 

32 

1 

1 

1  .026316  I 

.320000 

1  .  046715 

33 

1 

1 

1  .026316  I 

.320000 

1  .  046715 

34 

1 

1 

1  .026316  1 

.000000 

1  .OOOOOO 

36 

1 

4 

I  .105263  | 

.350000 

1 . 051095 

37 

1 

1 

1  .026316  | 

.350000 

1 .051095 

38 

1 

2 

I  .052632  | 

.350000 

1  .051095 

39 

1 

1  .026316  I 

.320000 

1 . 046715 

40 

1 

2 

1  .052632  | 

.  35C000 

1 . 051095 

41 

1 

1 

1  .026316  1 

.  400000 

1  .  058394 

43 

i 

1 

1  .026316  I 

. 350000 

1 . 051095 

45 

1 

1 

1  .026316  I 

. 350000 

1 . 051095 

46 

1 

1 

1  .026316  I 

. 45000C 

1  .065693 

48 

1 

1 

I  .026316  1 

. 350C00 

! . 051095 

56 

1 

1 

1  .026316  I 

. 300000 

!  . 043796 

62 

1 

X 

1  .026316  | 

. 3500C0 

1 . 051CS5 

68 

1 

1 

1  .026316  I 

. 53C0C0 

1  .  077372 

86 

1 

1 

1  .026316  | 

.50CC30 

!  .  072993 

70 


PMF 


cost 


.021  .042  .063  .084  .105 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 
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cost 


Membership  Function 
.106  .212  .318  .424  .530 


10  | 
11  I 

14  I 

15  I 
17  | 
20  | 
21  I 

23  | 

24  I 
27  | 

31  I 

32  ! 

33  | 

34  | 

36  | 

37  | 

38  1 

39  | 

40  | 

41  I 
43  | 

45  | 

46  | 
48  | 
56  | 
62  | 
68  1 
86  | 


cost 


Fuz2y  Probabilities 
.015  .031  .046  .062  .077 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


73 


APPENDIX  C 

Sensitivity  Analysis  Output  Data 


7-1 


OUTPUT  CORRESPONDING  TO  GOOD  LOCALITY 


THE  FINAL  RESULTS 


COST  I  No.  OF  OBS.  I  PMF  I  MEMBERSHIP  VALUES  I  FUZZY  PROBABILITIES 


10  | 

1 

1  .026316  I 

. 000000 

1 . 000000 

11  1 

1 

1  .026316  | 

. 000000 

1  .  000000 

14  | 

1 

1  .026316  I 

. 150000 

1  . 02051  * 

15  I 

2 

1  .052632  I 

. 150000 

1 . 020216 

17  1 

1  .026316  I 

. 000000 

1  .  000000 

20  I 

wBm 

wmm 

. 000000 

1  .000000 

21  1 

2 

1  .052632  | 

. 320000 

1  .  043127 

23  I 

3 

1  .078947  | 

. 600000 

1  .  080863 

24  I 

2 

1  .052632  | 

. 000000 

1 . 000000 

27  | 

1 

1  .026316  I 

.000000 

1  .  000000 

31  1 

1 

1  .026316  I 

.500000 

1  .067385 

32  1 

1 

1  .026316  I 

.450000 

1  .060647 

33  1 

1 

1  .026316  1 

.450000 

1  .060647 

34  1 

1 

1  .026316  I 

.450000 

1 . 060647 

36  I 

4 

1  .105263  1 

.350000 

1 .047170 

37  | 

1 

1  .026316  t 

.550000 

I . 074124 

38  | 

2 

1  .052632  I 

.320000 

1  .043127 

39  1 

1 

1  .026316  1 

.350000 

I  .047170 

40  I 

2 

1  .052632  I 

. 450000 

1  .060647 

41  1 

1 

1  .026316  1 

.320000 

I  .  043127 

43  I 

1 

1  .026316  I 

. 300000 

1 . 040431 

45  1 

1 

1  .026316  I 

.  320000 

I  .  043127 

46  | 

1 

1  .026316  1 

. 320000 

1 . 043127 

48  I 

1 

1  .026316  1 

.  150000 

1 . 020216 

56  i 

i 

1  .026316 

. 320000 

I  .  043127 

62  | 

1 

1  .026316  I 

. 150000 

1  .  020216 

68  1 

1 

1  .026316  l 

. 300000 

I . 04C431 

86  | 

1 

1  .026316  1 

.  150000 

I  .  020216 

75 


cost 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


PMF 

.021  .042  .063  .084  .105 

--I . I . I . I . I- 
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cost 


Membership  Function 
.120  .240  .360  .480  .600 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


Fuzzy  Probabilities 


cost 


•016  .032  .049  .065  . 0 e 1 


11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 
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OUTPUT  CORRESPONDING  TO  ION  PRICE 


THE  FINAL  RESULTS 


COST 

1 

No.  OF  OBS 

.  |  PMF  | 

MEMBERSHIP 

VALUES  I  FUZZ*  PROBABILITIES 

10 

1 

1 

I  .026316  | 

.000000 

1  .000000 

11 

1 

1 

I  .026316  | 

.000000 

1  .000000 

14 

1 

1 

I  .026316  I 

.150000 

1  .021368 

15 

1 

2 

1  .052632  | 

.150000 

1  .021368 

17 

1 

1 

1  .026316  | 

.100000 

1 .014245 

20 

1 

1 

I  .026316  | 

.000000 

1 .000000 

21 

1 

2 

I  .052632  I 

.670000 

I  .095442 

23 

1 

3 

1  .078947  | 

.670000 

1 .095442 

24 

1 

2 

I  .052632  | 

.000000 

I .000000 

27 

1 

1 

1  .026316  I 

.000000 

1 .000000 

31 

1 

1 

1  .026316  I 

.650000 

1 .092593 

32 

1 

1 

1  .026316  | 

.300000 

1.042733 

33 

1 

1 

1  .026316  I 

.450000 

1 .064103 

34 

1 

1 

1  .026316  I 

.450000 

I  .064103 

36 

1 

4 

1  .105263  I 

.350000 

1  .049858 

37 

1 

1 

I  .026316  I 

.3000O0 

I .042735 

36 

1 

2 

I  .052632  1 

.320000 

1 .045584 

39 

1 

8SSHP 

1  .026316  I 

.350000 

I  .049858 

40 

1 

2 

1  .052632  I 

.450000 

I  .064103 

41 

1 

1 

1  .026316  I 

.320000 

1  .  045584 

43 

1 

1 

1  .026316  | 

.150000 

i  .021368 

45 

1 

1 

I  .026316  1 

.300000 

1 .042735 

46 

1 

1 

I  .026316  I 

.300000 

1  .042735 

46 

1 

1 

1  .026316  | 

.040000 

I  .005698 

56 

1 

1 

1  .026316  I 

. 320000 

I  .045584 

62 

1 

1 

I  .026316  I 

.  040000 

| . 005698 

66 

1 

1 

1  .026316  I 

. 150000 

I  .021368 

66 

1 

1 

1  .026316  1 

. 040000 

1  .005698 

'  9 


PMF 


cost 


.021  .042  .063  .064  .105 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 
36 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


80 


cost 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


Membership  Function 
.134  .268  .402  .536  .670 


31 


Fuzzy  Probabilities 


cost 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


.019  .038  .057  .076  .095 


82 


OUTPUT  CORRESPONDING  TO  SIMPLE  TYPE 


THE  FINAL 

RESULTS 

COST 

1  No.  OF  OSS.  | 

PMF  |  MEMBERSHIP  VALUES 

FUZZY  PROBABILITIES 

10  I 

1  I 

.026316  | 

.  000000 

00000 

11  1 

1  1 

.026316  | 

.000000 

00000 

14  | 

1  I 

.026316  | 

.500000 

83333 

15  I 

2  1 

052632  I 

. 150000 

25000 

n  i 

. .  1  I 

.q-26316  | 

.000000 

000000 

20  | 

1  I 

.026316  I 

.000000 

000000 

21  I 

2  I 

.052632  I 

.320000 

053333 

23  1 

3  | 

.078947  | 

.  600000 

100000 

24  | 

2  I 

.052632  I 

. 000000 

000000 

27  | 

1  1 

.026316  1 

.000000 

000000 

31  1 

1  I 

.026316  | 

.  60000C 

100000 

32  | 

1  i 

.026316  i 

.220000 

036667 

33  I 

1  1 

.026316  I 

.220000 

036667 

34 

1  1 

.026316  I 

.220000 

036667 

36 

4  1 

.105263  I 

.170000 

028333 

37 

1  1 

.026316  I 

.620000 

103333 

38 

2  1 

.052632  | 

.220000 

036667 

39 

1  t 

.026316  | 

.170000 

028333 

40 

2  1 

.052632  | 

.220000 

036667 

41 

1  1 

.026316  | 

.  17000C 

028333 

43 

1  1 

.026316  t 

.220000 

036667 

45 

1  I 

.026316  I 

.320000 

053333 

46 

1  I 

.026316  1 

. 17C00C 

028333 

48 

1  1 

.026316  | 

.150000 

025000 

56 

1  1 

.026316  | 

.220000 

036667 

62 

1  1 

.026316  | 

. 150000 

025000 

66 

.1  1 

.026316  | 

.220000 

036667 

86 

1  1 

.026316  I 

. 150000 

025000 

83 


PMF 


cost 


.021  .042  .063  .084  .105 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
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OUTPUT  CORRESPONDING  TO  LOW  QUALITY 


THE  FINAL  RESULTS 


COST  I  No.  OF  OBS 

.  I  PMF  | 

MEMBERSHIP 

VALUES! FUZZY  PROBABILITIES 

10  1 

1 

i  .026316  I 

. 150000 

1  .016094 

11  1 

1 

1  .026316  I 

. 400000 

1  .  042918 

14  1 

1 

1  .026316  I 

. 150000 

1 . 016094 

15  I 

2 

1  .052632  I 

. 600000 

1 .064378 

17  | 

1 

I  .026316  I 

. 000000 

1  .  000000 

20  I 

1 

I  .026316  ! 

. 600000 

1  .  064378 

21  I 

2 

1  .052632  I 

.450000 

1  . 048283 

23  I 

3 

1  .078947  | 

. 620000 

1  .  066524 

24  1 

2 

!  .052632  I 

. 350000 

1 . 037554 

27  | 

1 

I  .026316  I 

.350000 

1  .037554 

31  1 

1 

I  .026316  | 

.500000 

1  .053648 

32  I 

1 

1  .026316  1 

.450000 

1  .048283 

33  I 

1 

I  .026316  1 

.450000 

1  .048283 

34  I 

1 

I  .026316  1 

.450000 

1 . 048283 

36  1 

4 

1  .105263  l 

.350000 

l .037554 

37  | 

1 

1  .026316  1 

.500000 

1  .053648 

38  1 

2 

1  .052632  I 
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1  .032189 

39  1 

1 

1  .026316  I 

.250000 

1  .  026824 

40  1 

2 

1  .052632  I 

.450000 

1  .048283 

41  1 

1 

I  .026316  I 

.250000 

1  .026824 

43  1 

1 

1  .026316  1 
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1  .  026824 
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1 
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OUTPUT  CORRESPONDING  TO  LOW  TECHNOLOGY 


THE  FINAL  RESULTS 

COST  I  No.  OF  OBS 

.  1  PMF  | 

MEMBERSHIP 

VALUES  I  FUZZY  PROBABILITIES 

10  | 

1 

1  .026316  | 

. 000000 

1 . 000000 

11  1 

1 

1  .026316  1 

. 000000 

1  .  000000 

14  : 

1 

1  .026316  I 

. 150000 

1  .  019711 

15  I 

2 

1  .052632  I 

. 150000 

1  .  019711 

17  | 

1 

1  .026316  | 

.000000 

1  .  oooooc 

20  | 

1 

1  .026316  | 

. 000000 

1  .000000 
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2 

1  .052632  | 

.320000 
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23  1 

3 

1  .078947  | 
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1 . 081472 

24  | 

2 

1  .052632  I 
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27  | 

1 

i  .026316  ! 

.000000 

1  .  000000 
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1  .026316  1 
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1  .070844 
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1 

1  .026316  | 

.450000 

1 .059133 
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1  .026316  | 

.450000 

1  .059133 
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1 

1  .026316  | 

.450000 

1  .  059133 
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.350000 
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1 

I  .026316  I 
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38  | 

2 
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1  .042050 
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1  .026316  1 

.350000 

1  .045992 

40  | 

2 

I  .052632  | 

. 450000 

1  .059133 

41  1 
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1  .026316  | 
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1  .042050 

43  1 
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1  .026316  1 
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THE  FINAL  RESULTS 


COST  1  No.  OF  OBS 

.  1  PMF  1 

MEMBERSHIP 

VALUES  1  FUZZY  P? 1BABILITIES 

10  i 

1 
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1 . 000000 
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14  ! 
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21  I 
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1  .078947  | 

. 620000 

1 . 081472 
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1  .052632  | 

.  000000 
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27  | 
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1  .  000000 

31  1 
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32  I 
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.450000 

I  .059133 

33  I 
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1  .026316  | 

.450000 

1  .059133 

34  | 

1 

1  .026316  | 

.450000 

1  .059133 

36  1 

4 

1  .105263  1 

. 350000 

1  .045992 

37  | 

1 

I  .026316  | 

. 620000 

I . 081472 

38  I 

2 

1  .052632  | 

.320000 

1 .042050 

39  1 

1 

1  .026316  | 

. 350000 

I  .045992 

40  I 

2 

1  .052632  | 
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1  .  059133 

41  1 

1 
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. 320000 

1  .  042050 
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1 

1  .026316  i 

. 150000 
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I  .026316  | 
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62  I 

1 

1  .026316  | 

.150000 
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68  1 

1 

I  .026316  | 
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1  .039422 
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OUTPUT  CORRESPONDING  TO  HIGH  PRICE 


THE  FINAL  RESULTS 

COST  |  No.  OF  OBS 

.  I  PMF  I 

MEMBERSHIP 

VALUES  1  FUZZY  PROBABILITIES 

10  | 

1 

1  .026316  I 

. 000000 

i . coco:: 

11  1 

1 

1  .026316  I 

. 000000 

1  .000000 

14  | 

1 

1  .026316  I 

.150000 

1  .  019763 

15  | 

2 

1  .052632  1 

. 150000 

1  .019763 

17  | 

1 

1  .026316  i 

.000000 

1 . 000000 

20  I 

1 

i  .026316  I 

.  000000 

1  .  000000 

21  1 

2 

1  .052632  1 

. 150000 

1 .019763 

23  I 

3 

1  .078947  | 

. 400000 

1  .  052701 

24  I 

2 

1  .052632  l 

. 000000 

1  .  000000 

27  | 

1 

1  .026316  I 

. 000000 

1  .  000000 

31  1 

1 

1  .026316  I 

.300000 

I  .039526 

32  I 

1 

1  .026316  I 

.450000 

1  .059289 

33  I 

1 

1  .026316  | 

.450000 

1  .  059289 

34  | 

1 

1  .026316  I 

.400000 

1  .052701 

36  I 

4 

1  .105263  I 

.350000 

I  .046113 

37  | 

1 

\  .026316  1 

.670000 

1  .088274 

38  I 

2 

1  .052632  I 

. 450000 

1  .059289 

39  1 

1 

I  .026316  1 

.350000 

1  .046113 

40  1 

2 

1  .052632  | 

. 450000 

1  .  059289 

41  1 

1 

1  .026316  1 

. 320000 

1 . 042161 

43  1 

1 

1  .026316  1 

. 320000 

1  .  042161 

45  1 

1 

1  .026316  I 

. 320000 

I  .  042161 

46  | 

1 

I  .026316  I 

. 320000 

1  .  042161 

48  | 

1 

1  .026316  I 

.  450000 

1  .  059289 

56  | 

1 

1  .026316  1 

.  300000 

1  .039526 

62  1 

1 

1  .026316  1 

.  320000 

1 . 042161 

68  I 

1 

1  .026316  I 

.  320000 

1  . 042161 

86  I 

1 

1  .026316  1 

.200000 

I  .  026350 
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OUTPUT  CORRESPONDING  TO  COMPLICATED  TYPC 


THE  FINAL  RESULTS 
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OUTPUT  CORRESPONDING  TO  HIGH  QUALITY 


THE  F 

INAL  RESULTS 

COST  1  No.  OF  OBS 

.  1  PMF 

.  MEMBERSHIP 

VALUES! FUZZY  PROBABILITIES 

10  1 

1 

1  .026316  1 

.  000000 

I  ri  ft  n  r>  n  ft 

1  -  w  J  ^  v  J 

11  1 

1 

1  .026316  1 

.  C  C  C  0  0  0 

!  .030000 

14  1 

1 

1  .026316  1 

.  150000 

i  .219894 

15  1 

2 

1  .052632  1 

.15  0 0  0 0 

!  .2  198  94 

17  | 

1 

!  .026316  ! 

•  V  w  «  W  J  J 

i  r-  f\  f\  *  ~ 

\  •  V  W  ^  w 

20  I 

1 

1  .026316  I 

.  000000 

I  .  0  0  0  0  0  c 

21  1 

2 

1  .052632  I 

. 32C000 

i  .042440 

23  1 

3 

i  .078947  ! 

. 320000 

I  .  042440 

24  | 

2 

1  .052632  1 

. 000000 

1  .000000 

27  | 

1 

1  .026316  1 

.000000 

1 . 000000 

31  1 

1 

1  .026316  1 

. 600000 

1  .  079576 

32  | 

1 

1  .026316  | 

. 320000 

1  .042440 

33  1 

1 

1  .026316  1 

. 32C000 

i  .  042440 

34  I 

1 

i  .026316  i 

.450000 

1  .059682 

36  i 

4 

i  .105263  i 

. 350000 

1  .046419 

37  | 

1 

1  .026316  1 

.  62 0 C 0 0 

1  032228 

38  I 

2 

1  .052632  1 

. 350000 

1  .046419 

39  i 

1 

1  .026316  1 

. 322000 

1  .  042440 

40  1 

2 

i  .052632  i 

. 450000 

i  .  059682 

41  1 

1 

!  .025316 

.350200 

346419 

43  1 

1 

a 

1  .026316  1 

.  v  0  0  0  c 

.239*88 

45  1 

1 

1  .026316  i 

. 62 00  CO 

!  .  082228 

46  1 
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1  .026316  1 

. 350000 
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43  | 
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. 150000 

. 019894 
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OUTPUT  CORRESPONDING  TO  HIGH  TECHNOLOGY 


THE 

FINAL  RESULTS 

COST  I  No.  or  OBS 

.  I  PMF 

1  MEMBERSHIP  VALUES 

IFU22Y  PROBABILITIES 

10  I 

1 

I  .026316 

1  .000000  1 

.  000000 

11  1 

1 

I  .026316 

1  .000000  1 

.  000000 

14  I 

1 

1  .026316 

1  .150000  1 

. 0233  64 

15  1 
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n mm 
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. C23364 

17  | 

1  .026316 

1  .000000  1 

.000000 

20  I 

■ 

1  .000000  1 

.  000000 

21  1 

2 

1  .052632 

1  .320000  1 
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23  I 

3 

1  .078947 

1  .350000  1 

.  054517 

24  I 

2 

I  .052632 

1  .000000  1 

.  000000 

27  | 

i  .026316 

I  .000000  1 
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|  .026316 
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33  I 

i 

I  .026316 

1  .350000  I 

.  054517 

34  I 

I  .026316 

1  .350000  1 

.  054517 

36  | 

1  .105263 

!  .350000  ; 

.  054517 

37  | 

HI 

1  .026316 
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1 

1  .026316 

!  .350000  1 

.054517 

40  1 

2 
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1  .320000  ! 

. 049844 

48  1 

1 

I  .026316 

1  .150000  I 

. 023364 

56  1 

1 

1  .026316 

!  .320000  1 

. 049844 

62  1 

1 

1  .026316 

l  .150000  ! 

. 023364 

63  I 

1 

1  .026316 

i  .30CC00  i 

.  C  4  6"?  2  9 

86  I 


.  C26316 


.  150000 


!  .  02336  . 


1  l  i. 


Membership  Function 


cost 


10 

11 

14 

15 
1? 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 


.070  .140  .210  .280  .350 

"I - I . I - I - I- 


1  i  3 


Fuzzy  Probabilities 


cost 


.011  .022  .033  .044  .055 


10 

11 

14 

15 
17 
20 
21 

23 

24 
27 

31 

32 

33 

34 

36 

37 

38 

39 

40 

41 
43 

45 

46 
48 
56 
62 
68 
86 
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APPENDIX  D 


User  Manual 


User  Manual 


1  Introduction 

The  preliminary  cost  estimate  is  a  primary  consideration  of  both  owners  and 
designers  of  construction  projects.  The  estimation  at  this  stage  involves  many 
design  variables.  Usually  the  designer  does  not  have  enough  engineering  data 
to  predict  the  future  behavior  of  these  variables  other  than  verbal  description, 
such  as  ’  good.’’  "bad.'  "large.’’  ’  small.’  3nu  so  on. 

Using  regression  analysis,  this  program  uses  previous  data  from  the  past 
projects  to  develop  a  cost  function. 

The  opinion  of  the  experienced  designer  regarding  the  value  of  the  design 
variables  is  u^ed  to  find  cost  membership  function. 

The  variable  membership  values  range  from  0  to  1  depending  upon  the 
degree  of  confidence  that  the  variable  is  a  member  of  the  set.  Applying 
the  extension  principle,  membership  functions  of  independent  design  vari¬ 
ables  are  used  to  find  the  cost  membership  function.  The  computer  program 
(COSTEST)  developed  in  this  study  enables  the  estimator  to  obtain: 

1  A  membership  function  which  explains  the  degree  that  a  cost  value 
belongs  to  a  c:>st  event. 

2.  The  probability  of  each  cost  value  and  event.  (The  second  one  is  com¬ 
puted  using  the  probability  mass  function  (PMF)  of  the  cost.) 

An  example  is  given  below  to  show  the  input  and  output  of  the  computer 
program. 


2  Example  Problem 

Assume  that  the  independent  design  variables  are  \  \ .  U2.  1  3.  and  1  \  and  C  is 
the  dependent  design  variable  representing  the  unit  cost  of  the  project.  The 
linear  function: 

C  —  C2o  —  ~  0-2^2  ~  a3^3  —  ^-1^4  (1) 

represents  the  relation  between  the  independent  and  dependent  design  van- 


Table  1:  Historical  Data 


c, 

Vn 

r2l 

1  3: 

u. 

55.50 

4.5 

cc 

10 

*)  *T 

5.2 

58.00 

5.2 

4.2 

3.0 

4  0 

54.00 

4.7 

3.2 

1.5 

6.5 

64.00 

5.8 

4.9 

1.6 

6.9 

60.00 

4.0 

4.6 

2.5 

5.7 

49.00 

4.9 

3.0 

1.9 

4.8 

57.00 

4.5 

3.9 

2.8 

4.7 

57.00 

5.7 

4.2 

1.3 

6.3 

61.00 

5.3 

4.6 

2.4 

5  4 

54.00 

4.1 

3.2 

2.7 

5.0 

Table  I  shows  the  values  of  dependent  design  variable  C,  and  independent 

design  variables  \'kl.  k  =  (1 . 4).  i  =  (1,...,10)  associated  with  ten  past 

projects 

These  variables  can  be  indices  or  unit  prices  of  any  component  per  unit 
of  the  project.  This  unit  can  be  in  cubic  feet,  square  feet  or  other  units,  such 
as  cost  per  output  product,  for  example,  cost  per  student  in  school  projects. 

The  next  step  is  to  develop  a  membership  function  for  each  design  vari¬ 
able.  This  will  be  based  on  the  experienced  designer  opinion  of  the  behavior 
of  these  design  variables  in  the  future.  These  membership  functions  can  be 
discrete  or  continuous.  In  this  example,  assume  the  following  membership 
functions  for  the  four  design  variables. 

Assume  \\  represents  the  weather  index.  represents  the  productivity 
index.  \  3  the  price  index  and  \\  the  quality  index. 

Suppose,  in  this  example,  the  prediction  of  design  variables  is  as  shown 
in  Table  2. 

Figures  1,  2,  3  and  4  represent  the  assumed  membership  functions  of 
variables  V1(  \ \ ,  1  3  and  \’4,  respectively. 

According  to  the  designer  opinion  in  Table  2  for  the  prediction  of  design 
variables  \  \  (see  Table  2)  and  using  Figures  1-4  and  Table  1.  the  membership 
values  of  these  design  variables  are  computed  and  shown  in  Table  3.  How  to 
use  these  membership  values  is  explainer!  m  the  following  section. 
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mtxm.  moo 


Figure  2:  Productivity  Index  Membership  Function 


Figure  3:  Price  Index  Membership  Function 


Figure  4  Quality  Index  Membership  Function 
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Table  2:  Designer  Expectation  of  the  Variables 


Variable 

v; 

Future  Prediction  . 

Weather  Index 

V, 

Good  1 

Productivity  Index 

Vj 

Medium 

Price  Index 

r3 

High  1 

Quality  Index 

High  : 

3  Input 

This  program  is  developed  to  read  the  data  from  a  file.  The  program  indicates 
the  name  of  the  input  file. 


3.1  Input  File 

The  input  file  should  contain  the  dependent  design  variable  values,  C,,  in¬ 
dependent  design  variable  values.  I*,-,  and  membership  values.  in  the 
format  as  shown  in  Table  4. 

Once  the  input  file  is  established,  the  program  can  be  run.  The  program 
contains  an  interactive  part  to  prompt  the  user  to  input  some  information, 
as  in  the  following  steps,  while  running  the  example  problem. 

■  a  >  costest; 

To  run  the  program  from  the  PC. 


'  SRun  COSTEST  j 

To  run  the  program  from  the  VAX. 


,  Please  enter  the  output  title  max  40  characters. ' 

Prompts  the  user  to  input  the  output  files  title.  This  title  will  appear  at 
the  top  of  all  output  files, 


Example  Problem 

This  is  the  title  entered  by  the  user. 


Please  enter  the  input  file. 
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Table  3:  Design  Variables  and  Their  Corresponding  Membership  Values 


:  C‘ 

Vu 

0  v„ 

( \  i ,  =  Good) 

l  2i 

Mv2, 

(l'2l  =  Med) 

^  3i 

Os,  =  High) 

U, 

(U,  =  High) 

55.50 

4.5 

.75 

3.2 

2 

2.2 

0.6 

5.2 

.433 

:  58.00 

5.2 

.4 

4.2 

.8 

3.0 

1.0 

4.0 

0.0 

54.00 

4.7 

.65 

3.2 

2 

1.5 

.25 

6.5 

.833 

64.00 

5.8 

.1 

4.9 

.1 

1.6 

.3 

6.9 

.967 

:  60.00 

4.0 

1.0 

4.6 

.4 

2.5 

.75 

5.7 

.567  1 

49.00 

4.9 

.55 

3.0 

0.0 

1.9 

.45 

4.8 

.267 

57.00 

4  5 

.75 

3.9 

.9 

2.8 

.9 

4.7 

.233  ; 

57.00 

5.7 

.15 

4.2 

.8 

1.3 

.15 

6.3 

.767 

61.00 

5.3 

.35 

4.6 

.4 

2.4 

.  1 

5.4 

.467 

54.00 

4.1 

0.95 

3.2 

.20 

2.7 

.85 

5.0 

.333 

Prompts  the  user  to  enter  the  input  file  that  is  explained  in  the  input 
section  which  contains  the  dependent  design  variables,  independent  design 
variables  and  membership  function  values. 


EXCOST.DAT 

This  is  the  name  of  the  input  file. 


Please  enter  regression  analysis  output  file 

Prompt  the  user  to  enter  output  file  that  contains  the  regression  matrix, 
cost  function  coefficient,  quantification  of  error  of  linear  regression  and  the 
final  cost  function. 


REGEX  OUT 

This  file  is  the  regression  output  file  for  the  example. 

Please  enter  probabilities  and  membership  functions  output  file 
Prompts  the  user  to  enter  the  output  file  that  contains  cost  function, 
table  that  contains  cost,  number  of  observations  of  that  cost,  probability  mass 
function,  and  corresponding  membership  function  and  cost  fuzzy  probability. 
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This  program  will  generate  three  plots  They  are  probability  mass  func¬ 
tion  (PMF).  cost  membership  function  and  cost  fuzzy  probability.  All  of 
three  plots  will  also  be  in  this  output  file. 


CMEX.Ol'T 

This  is  the  name  of  the  example  output  entered  by  the  user. 


Please  enter  data  output  file 

Prompts  the  user  to  input  an  output  file,  which  contains  the  input  data, 
for  review. 


IN  EX  OUT 
This  is  the  output  file. 


Plot  option. 


Pleas  e  enter  (1)  if  plot  is  needed. 

(2)  if  plot  is  not  needed. 

As  mentioned  earlier,  this  program  will  generate  three  plots,  but  these 
plots  are  optional.  So  the  program  will  prompt  the  user  to  choose  between 
plotting  or  not.  If  (2)  is  selected,  the  program  will  stop.  But  if  (1 )  is  selected, 
then  the  following  message  will  appear: 


Line  max  length  option 

Please  enter  (1)  to  use  default  (max.  L  =  40) 

’  (2)  to  use  new  value  for  max.  L 

This  option  is  to  chose  the  maximum  line  length  in  the  plot.  If  (1)  is 
selected,  the  program  will  stop  and  the  maximum  line  in  the  plot  will  be 
40  and  the  other  lines  will  be  scaled  accordingly  But  if  (2)  is  chosen,  the 
program  will  ask  for  a  value. 


Please  enter  line  max  length 
°30~ 

30  will  be  the  maximum  length  in  the  plots  and  as  in  (1)  the  other  values 
will  be  scaled  accordingly 
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Table  4:  Inpt  Data  File  Format 


A- 

1- 

Number 

of  Variables 

B- 

1- 

Number 

of  Data  Points 

C- 

1- 

Ci  r„ 

1  12 

1  13  . 

^  1  m 

o_ 

^ii 

Ml2 

A*  13  . 

AM  m 

3- 

C2  l*2i 

\  22 

c  23  . 

^  2m 

4- 

P21 

^  23  . 

H  2m 

2n  —  1-  Cn  \  ni 

1’nJ 

l'n3  ••• 

^  nm 

2  n-  Mm 

Mn2 

fJ-n3  ■■■ 

A*  nm 

The  three  computer  output  files  are  shown  in  the  next  section. 


4  Output 

The  output  is  contained  in  three  files,  as  explained  in  the  previous  section. 
Computer  output  of  these  three  files  is  added  to  this  manual  for  illustration. 
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OUTPUT  FILES 


Rgression  Analysis  Output  File 
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EXAMPLE  PROBLEM 


OUTPUT 


10 . ooooco 

48 .700001 

39 . 000000 

21 . 900000 

54.50000 

4 

5S9. 500000 
48.700001 

240 . 669998 

191 .790009 

104 . 940002 

267 .27002 

: 

2 '’32 . 549805 
39.000000 

191 .790009 

156.540009 

85 . 540009 

214 .30000 

3 

2045 . 700195 
21.900000 

104 . 940002 

85 . 540009 

51.090000 

115.41000 

4 

1248 .500000 

54 . 50C004 

267.270020 

214 .300003 

115 .410004 

304.3-002 

6 

mm 

3116 . 80CC49 

mmm mm  *  w  raw  mmmmm 

mm mm^ mm m m m m m m m 

1 

1 

i 

1 

l 

1 

l 

s 

1 

1 

1 

1 

• 

1 

COST  FUNCTION  COEFFICIENTS 


a  0-  -3.817445 
a  1-  2 . 3307 9.5 
a  2-  2.661444 
a  3-  7.042204 
a  4-  4.332922 


QUANT IF I CAT ION  OF  ERROR  OF  LINEAR  REGRESSION 


THE 

THE 

THE 

THE 

THE 


TOTAL  Or  SUM  OF  SQUARES  Or  RESIDUALS  St-  32592.250 


SUM  OOF  SQUARES  Or  TH  RESIDUALS 
STANDEE  ERROR  OF  ESTIMATE 
COEFFICIENT  OF  DETERMINATION 
CORRELATION  COEFFICIENT 


Sr-  6.889 

Sx/y-  1.174 

r'2-0 . 9998 
r—0 . 9999 


THE  FINAL  COST  FUNCTION 

y-  -3 . 81745*  ( 

2.33079)  X  l+(  2.66144)  X  2+  (  7.04220)  X  3* (  4  332  92) 
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Membership  Function  and  Fuzzy  Probability 
Output  File 
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Membership  Function  and  Fu2zy  Probability 
Output  File 
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THE  COST  FUNCTION 


EXAMPLE  PROBLEM 


vm  -3 . 81745+  (  2.33079)  X  1+ (  2.66144)  X  2*  (  7.04220)  X  3+ (  4.33292 
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THE  COSTS  AND  THEIR  CORRESPONDING  MEMBERSHIP  FUNCTIONS 


COST 

|  MEMBERSHIP  VALUES 

1 

50 

I  0.00000 

1 

53 

I  0.20000 

i 

54 

1  0.20000 

1 

55 

I  0.20000 

1 

57 

1  0.23300 

1 

57 

I  0.15000 

1 

58 

1  0.00000 

1 

60 

|  0.40000 

1 

61 

I  0.35000 

1 

64 

!  0.10000 

1 
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THE  riNAL  RESULTS 


COST 

1 

Ho.  OF  OBS. 

PMT 

1 

MEMBERSHIP  VALUES 

50 

1 

1 

.100000 

1 

.000000 

53 

1 

1 

.100000 

1 

.200000 

54 

1 

1 

.  100000 

1 

.200000 

55 

1 

1 

.  100000 

1 

.200000 

57 

1 

2 

.200000 

1 

.233000 

58 

1 

1 

.100000 

1 

.000000 

60 

1 

1 

.  100000 

1 

.400000 

61 

1 

1 

.100000 

1 

.350000 

64 

1 

1 

.  100000 

1 

.100000 

50 

53 

54 

55 

57 

58 
60 
61 
64 


Membership  Function 

.080  .160  .240  .320  .400 

—  I . I . I . I . I- 
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cost 


50 

53 

54 

55 

57 

58 
€0 
Cl 
C4 


russy  frsb*biliti«« 

.048  .095  .143  .190  .238 

— , - 1 - 1 - 1 - I- 
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